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ABSTRACT 


We characterize the earliest galaxy population in the JADES Origins Field (JOF), the deepest 
imaging field observed with JWST. We make use of the ancillary Hubble optical images (5 filters 
spanning 0.4—0.9um) and novel JWST images with 14 filters spanning 0.8 — 5um, including 7 medium- 
band filters, and reaching total exposure times of up to 46 hours per filter. We combine all our data 
at > 2.3um to construct an ultradeep image, reaching as deep as ~ 31.4 AB mag in the stack and 
30.3-31.0 AB mag (50, r = 0.1” circular aperture) in individual filters. We measure photometric 
redshifts and use robust selection criteria to identify a sample of eight galaxy candidates at redshifts 
z = 11.5 — 15. These objects show compact half-light radii of Ry/2 ~ 50 — 200pc, stellar masses of 
M, ~ 107 — 10* Mo, and star-formation rates of SFR ~ 0.1—1 Mo yr !. Our search finds no candidates 
at 15 « z « 20, placing upper limits at these redshifts. We develop a forward modeling approach to 
infer the properties of the evolving luminosity function without binning in redshift or luminosity that 
marginalizes over the photometric redshift uncertainty of our candidate galaxies and incorporates the 
impact of non-detections. We find a z — 12 luminosity function in good agreement with prior results, 
and that the luminosity function normalization and UV luminosity density decline by a factor of ~ 2.5 
from z = 12 to z = 14. We discuss the possible implications of our results in the context of theoretical 
models for evolution of the dark matter halo mass function. 
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1. INTRODUCTION 


JWST has pushed the forefront of our knowledge of 
galaxies in the distant universe to the first 350 million 
years of cosmic time. Within the first weeks of oper- 
ations, surveys with JWST unveiled galaxy candidates 
beyond redshift z ~ 12 in an epoch when only the most 
optimistic models of the cosmic star formation rate den- 
sity predicted that galaxies would be easily discoverable 
(Naidu et al. 2022a; Castellano et al. 2022; Finkelstein 
et al. 2023a; Adams et al. 2023a; Atek et al. 2023; Don- 
nan et al. 2023b; Harikane et al. 2023b; Morishita & 
Stiavelli 2023; Bouwens et al. 2023). The identifica- 
tion and spectroscopic confirmation by the JWST Ad- 
vanced Deep Extragalactic Survey (JADES; PIs Rieke 
and Lutzgendorf; Eisenstein et al. 2023a) of the galaxies 
JADES-GS-z12-0 at z = 12.6 and JADES-GS-z13-0 at 
z = 13.2 affirmatively established for the first time the 
reality of galaxies at z > 12 (Curtis-Lake et al. 2023; 
Robertson et al. 2023; D'Eugenio et al. 2023). Subse- 
quently, other galaxy candidates have been confirmed at 
z ~ 12—13 in other surveys (Fujimoto et al. 2023; Wang 
et al. 2023a) and many additional high-redshift can- 
didates identified photometrically (e.g., Hainline et al. 
2023a; Pérez-González et al. 2023a; Leung et al. 2023). 

The discovery of these distant sources raises sub- 
stantial questions about the nature of galaxy forma- 
tion in the early universe (Ferrara et al. 2023; Mason 
et al. 2023; Dekel et al. 2023; Li et al. 2023; Lovell 
et al. 2023; Shen et al. 2023a; Yung et al. 2024). The 
earliest known galaxies appear relatively bright (e.g., 
Naidu et al. 2022a; Castellano et al. 2022; Treu et al. 
2023; Finkelstein et al. 2023b), show a range of stel- 
lar masses M, ~ 10° — 10? Ma, and have young stel- 
lar ages of tą ~ 10/ — 108 yr (Robertson et al. 2023). 
Structurally, these galaxies show physical sizes of r ~ 
0.1 — 1 kpc and star formation rate surface densities of 
X, 50 — 100 Mo yr-lkpc ? (Robertson et al. 2023; 
Arrabal Haro et al. 2023; Wang et al. 2023a). They are 
compact star forming galaxies undergoing rapid star for- 
mation on a timescale comparable to their local dynam- 
ical times. Individually, the properties of these objects 
are not extreme given the densities and dynamics of the 
early universe. Collectively, the apparent, albeit uncer- 
tain, abundance of such objects in the context of struc- 
ture formation may be unexpectedly high. Resolving 
this essential quandary requires statistical constraints 
on the abundance of z > 12 galaxies and information on 
their possible origins through higher-redshift searches. 


To answer these questions, this work presents first re- 
sults on the search for distant galaxies in the JADES 
Origin Field (JOF; Program ID 3215, PIs Eisenstein 
and Maiolino; Eisenstein et al. 2023b). The JOF obser- 
vations were designed to use JWST medium bands, in- 
cluding NIRCam F162M, to isolate the Lyman-o break 
at z = 12 and simultaneously control for contamina- 
tion by lower-redshift line emitters that can mimic the 
broad-band spectral energy distributions (SEDs) of dis- 
tant galaxies (Naidu et al. 2022b; Zavala et al. 2023; 
Arrabal Haro et al. 2023; Pérez-González et al. 2023b). 
In concert with ultra-deep broad-band observations from 
JADES, the 9.05 arcmin? JOF provides the best current 
dataset for finding and characterizing z = 12 galaxies. 
We search the JOF for objects to an effective limiting 
depth of f, ~ 2 — 3 nJy, performing SED fitting anal- 
yses to select the highest redshift candidates. We then 
use a forward-modeling approach to infer the character 
of the evolving luminosity function given the properties 
of our sample of high-redshift candidate galaxies. Our 
method accounts for the photometric redshift posterior 
constraints of our sample's galaxies without binning in 
redshift or luminosity. We employ our method to study 
the behavior of the evolving luminosity function beyond 
z ~ 12 and the abundance of galaxies at earlier times. 

This paper is organized as follows. In 82 we review the 
JOF data, the observations, data reduction procedure, 
source detection, and photometry. In 83, we describe 
our selection procedure based on SED template fitting. 
Forward modeling constraints on the galaxy candidate 
structural properties and inference of the distant stellar 
population properties are described in 84. We charac- 
terize the galaxy luminosity functions at z ~ 12 — 15 
and our constraints on the UV luminosity density at 
z ~ 12 — 20 in 85, and report the inferred physical prop- 
erties of the high-redshift candidates in 86. We interpret 
the observational results in the context of galaxy forma- 
tion theory in 87. We summarize our conclusions and 
preview future work in 88. Throughout this work, we 
use the AB magnitude system (Oke & Gunn 1983) and 
assume a flat Lambda cold dark matter (ACDM) cos- 
mology with Qm = 0.3 and Ho = 70 km s^! Mpc !. 


2. DATA 


This work uses JWST observations in the JOF to dis- 
cover and constrain the abundance and properties of 
z > 12 galaxies. In 82.1 we review the JOF and ac- 
companying JADES and Hubble Space Telescope (HST) 


observations. In 82.2, we present the data reduction 
methods used to process the imagery. The detection 
and photometric methods used to discover the objects 
are described in 82.3. 


2.1. Observations 


Eisenstein et al. (2023b) presents the JOF, a single 
JWST NIRCam pointing of exceptional depth, with 
about 7 days of exposure time spread between 14 fil- 
ters covering an A ~ 9 arcmin? area. The JOE began 
with the parallel imaging of deep JADES spectroscopy 
(Program ID 1210, presented in Bunker et al. 2023) 
that produced long F090W, F115W, F150W, F200W, 
F277W, F335M, F356W, FA10M, and F444W exposures 
in a field adjacent to the Hubble Ultra Deep Field within 
the GOODS-S field. This campaign continued in Cycle 
2 Program ID 3215, which observed in 6 JWST NIR- 
Cam medium bands—F162M, F182M, F210M, F250M, 
F300M, and F335M—again acquired in parallel to deep 
NIRSpec observations. We also include all JADES 
GOODS-S medium-depth imaging (Program ID 1180) 
that overlaps with the JOF. This area of GOODS-S 
partially overlaps with the FRESCO (Program ID 1895) 
F182M, F210M, and F444W data, which we incorporate. 
The field also has partial coverage of HST ACS F435W, 
F606W, F775W, F814W and F850LP images reduced 
and released through the Hubble Legacy Field program 
(Illingworth et al. 2016a) reductions of the Great Obser- 
vatories Origins Deep Survey (GOODS; Giavalisco et al. 
2004) and Cosmic Assembly Near-infrared Deep Extra- 
galactic Legacy Survey (CANDELS; Grogin et al. 2011; 
Koekemoer et al. 2011) images. In total, these data pro- 
vide nineteen JWST and HST photometric bands that 
we use to constrain the galaxy SEDs and particularly 
the Lyman-a break. 


2.2. Data reduction 


Our image reduction methods were outlined in Rieke 
et al. (2023) and Eisenstein et al. (2023a), detailed in 
Tacchella et al., (in prep), and we provide a summary 
here. We process the images with the jwst Calibration 
Pipeline (version 1.11.4) and Calibration Reference Data 
System pipeline mapping (CRDS pmap) 1130, which 
includes in-flight NIRCam dark, distortion, bad pixel 
mask, read noise, superbias and flat reference files. 

We use jwst Stage 1 to perform the detector- 
level corrections and ramp fitting. We run this 
stage with the default parameters, except for the 
correction of “snowball” artifacts from cosmic rays. 
The identification and correction of snowballs rep- 
resent a big challenge. Heuristically, we find that 
the following parameters provide reasonable snow- 
ball amelioration: max jump.to flag neighbors = 
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1, min jump to flag neighbors = 100000, 
min_jump_area = 9, min sat area = 1, 
expand_factor = 2, min_sat_radius_extend = 2.5, 


and max-extended_radius = 200. 

As detailed in Rieke et al. (2023), we run jwst Stage 2 
with the default parameters, but replace the STScI flats 
for all LW bands except F250M and F300M with custom 
super-sky flats. When we do not have sufficient images 
to produce a robust flat field, we interpolated the flat- 
field images from the bands adjacent in wavelength. Fol- 
lowing Stage 2, we perform custom corrections for all ad- 
ditive effects including 1/f noise, scattered light effects 
(^wisps" and “claws”), and the large-scale background. 
Furthermore, we also updated the DQ data quality ar- 
ray to mask additional features imprinted visually onto 
the mosaics, including persistence, uncorrected wisp fea- 
tures, and unflagged hot pixels. 

Before running jwst Stage 3, we perform astrometric 
registration to Gaia DR2 (G. Brammer priv. comm., 
Gaia Collaboration et al. 2018) with a modified jwst- 
pipeline tweakreg code. We apply both a rotation 
and offset to the individual level-2 images. For images 
taken in the A module with the medium bands F182M, 
F210M, and F335M, we replace the default distortion 
maps with the nearest (in effective wavelength) wide- 
band distortion map for that detector. 

We construct the mosaics using jwst Stage 3. We cre- 
ate single mosaics for each filter by combining exposures 
from all observations, and run jwst Stage 3 with the de- 
fault parameter values while setting the pixel scale to 
0.03 ” /pixel and a drizzle parameter of pixfrac = 1 for 
the SW and LW images. Finally, we perform a cus- 
tom background subtraction, following the procedure 
outlined in (Bagley et al. 2023a). For F090W, F115W, 
and F150W, hot pixels that pass median rejection are 
replaced with median filtered values from the local flux 
image. 


2.9. Detection and Photometry 


The detection and photometry methods are discussed 
in Rieke et al. (2023) and Eisenstein et al. (2023a) and 
will be detailed in Robertson et al. (in prep). 

To perform source detection, an inverse variance- 
weighted stack of the long-wavelength NIRCam F250M, 
F277W, F300M, F335M, F356W, F410M, and F444W 
SCI and ERR channels are constructed. Small-scale 
noise residuals from incomplete masking in the must 
pipeline are median filtered from the ERR images. The 
signal-to-noise ratio (SNR) image created from the ratio 
of these images is used as the detection image. An initial 
source significance threshold of SNR > 1.5 is used to 
select regions of interest, and a series of custom compu- 
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tational morphology algorithms inspired by NoiseChisel 
(Akhlaghi & Ichikawa 2015; Akhlaghi 2019) are applied 
to refine the segmentations. Stars and diffraction spikes 
are masked by constructing segmentations from stacks of 
all available filters and integrated into the detection seg- 
mentation map. The detection image segmentations are 
deblended using a logarithmic scaling of the F200W im- 
age. High-pass filtering is applied to the outer regions of 
large segmentations to isolate proximate satellite galax- 
ies. After these refinements of the segmentation map, 
a final pass to detect potentially missed compact, faint 
sources is applied. The completeness as a function of 
flux and size for this detection algorithm has been calcu- 
lated using source injection simulations and is presented 
in Section 4.1. 

After engineering the segmentation map, we perform 
a set of customized photometric measurements based 
on the photutils (Bradley et al. 2023) analysis pack- 
age. Object centroids are computed using the “win- 
dowed positions" used by Source Extractor (Bertin & 
Arnouts 1996). Apertures for measuring Kron (1980) 
fluxes are determined based on the stacked signal im- 
age (the numerator of the SNR detection image) using 
a Kron parameter of 2.5. We limit the area of the 
Kron aperture to be less than twice an object's seg- 
mentation area. In addition to Kron fluxes, we mea- 
sure circular aperture photometry with aperture radii 
of r = {0.1”, 0.15”, 0.25”, 0.3”, 0.35”, 0.5” }. To provide 
aperture corrections, we produce a model point spread 
function (mPSF) following the method of Ji et al. (2023), 
where we inject WebbPSF models into jwst level-2 im- 
ages and mosaic them using the same exposure pat- 
tern as the JOF observations to provide a composite 
star field. An mPSF for each band (and observing pro- 
gram) is then constructed from these PSF-mosaics. The 
circular aperture corrections are measured and tabu- 
lated, and the Kron aperture corrections computed by 
integrating within the corresponding elliptical apertures 
placed on the mPSF. For HST, we measure empirical 
(e)PSFs using the photutils (Bradley et al. 2023) ePSF 
Builder with visually inspected stars in the field. 

We perform a bevy of photometric validation tests. 
Cross validation against the CANDELS survey HST 
photometry using completely independent HST re- 
ductions from the Hubble Legacy Field program are 
presented in Rieke et al. (2023) for the broader 
JADES/GOODS-S field. We also compute median 
photometric offsets from SED templates using EAZY 
(Brammer et al. 2008), following the method described 
by Hainline et al. (2023a). We find these zeropoint off- 
sets to be within 5.2%, and typically within 1%, for 
JWST filters. 


2.3.1. Surface Brightness Profile Modeling 


We also forward model each galaxy’s surface bright- 
ness profile using the Forcepho code (B. Johnson, in 
prep). We use Forcepho with custom model point- 
spread functions to model the surface brightness profile 
of each galaxy in our survey simultaneously with any 
nearby objects in each individual exposure where pixel 
covariance is minimized. We restrict the modeling to 
the F200W and F277W bands, to minimize the chance 
of any PSF mismatch or astrometric errors while maxi- 
mizing S/N and resolution. The surface brightness pro- 
file is assumed to be a Sérsic (1968) model, with a fast 
Gaussian-based factorization of the model. Forcepho 
provides a Bayesian estimate of the surface brightness 
profile parameters, including the galaxy half-light ra- 
dius. We have used Forcepho to study the structure 
of other extremely high-redshift galaxies (e.g., Robert- 
son et al. 2023; Tacchella et al. 2023), and we refer the 
reader to Baker et al. (2023) for more details on our 
morphological analysis methods. 


2.4. Image Depths 


With the construction of our broad- and medium-band 
NIRCam mosaics and the long-wavelength (A > 2.3m) 
detection image, we can use the photometry method de- 
scribed in §2.3 to measure our image depths. In Ta- 
ble 1, we report the median aperture corrected 5 — a 
point-source depth in each filter and the stack (using 
the F277W PSF to estimate the stack's aperture cor- 
rection). When measuring the depth in each image, we 
use a dilated version of the segmentation map created 
by the detection algorithm to mask source pixels. We 
note that the single-band images depths listed in Table 
1 are all within 10—25% of the 5—o point-source depths 
we reported in Eisenstein et al. (2023b) that were com- 
puted from the JWST exposure time calculator, with the 
longest wavelength filters showing the most improved 
depth. Our single-band images reach 30.3—31.0 AB, 
and the combined A > 2.3um stack reaches 31.4 AB 
depth. For comparison, we also list the 5 — o point- 
source depth the corresponding A > 2.3um stacks from 
available NIRCam long-wavelength images in NGDEEP 
(F277W+F356W+F444W; Bagley et al. 2023b), the 
MIRI-UDF NIRCam parallel (F277W+F356W; Pérez- 
González et al. 2023a), and the JADES GOODS-S Deep 
region that covers the Hubble Ultra Deep Field (Rieke 
et al. 2023). To measure their depths, we processed 
these fields using identical methods and used the same 
F277W PSF model to aperture correct them. We re- 
port depths for each program separately, and note that 
where the MIRI-UDF parallel and NGDEEP NIRCam 


Table 1. Depths of the JADES Origins Field 


Band Median Depth^ Median Depth 
[nJy] [AB] 
JWST/NIRCam Filters 
F090W 2.80 30.28 
F115W 2.99 30.48 
F150W 2.19 30.55 
F162M 2.76 30.30 
F182M 1.78 30.77 
F200W 2.21 30.51 
F210M 2.29 30.50 
F250M 2.58 30.37 
F277W 1.42 31.02 
F300M 1.80 30.76 
F335M 1.70 30.82 
F356W 1.58 30.90 
F410M 2.65 30.34 
F444W 2.26 30.52 
Stacked Depth at A > 2.3um 
JOF 0.96 31.44 
NGDEEP^* 0.82 31.61 
MIRI-UDF* 1.28 31.13 
JADES GOODS-S Deep 1.39 31.04 


NOTE—^ Median r = 0.1" aperture corrected 5e point-source 
depth. ^ This depth reflects our independent processing of 
the NGDEEP data, and we refer the reader to Bagley et al. 
(2023b) for their depth measurements. © This depth reflects 
our independent processing of the MIRI-UDF data, and we 
refer the reader to Pérez-González et al. (2023a) for their depth 
measurements. 


imaging overlap the combined depths will be even more 
sensitive than listed in Table 1. 


3. SELECTION OF REDSHIFT Z Z 12 GALAXIES 


The photometric selection of high-redshift galaxies re- 
lies on identifying a strong Lyman-o break in the rest- 
frame UV of a galaxy's SED (e.g., Guhathakurta et al. 
1990; Steidel et al. 1995). Below, we detail our selection 
of z = 12 galaxies based on this feature. 


3.1. Photometric Redshift Estimation 


To infer the photometric redshifts of galaxies in the 
JOF, we apply the techniques detailed in Hainline et al. 
(2023a) to fit templates of galaxy SEDs to our JWST 
and HST photometry, varying the redshift to assess the 
relative goodness of fit. To perform the SED fits, we 
use the EAZY code (Brammer et al. 2008) to compute 
rapidly the photometric redshift posterior distributions 
for each galaxy in the JOF survey. When fitting SED 
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templates, we use the template suite described in Hain- 
line et al. (2023a) that includes models with strong line 
emission and a range of UV continua. The photomet- 
ric redshifts estimated from fits to these templates were 
shown to have an outlier fraction (defined as the frac- 
tion of sources with |Zphot — Zspec|/(1 + Zspec) > 0.15) 
of fout = 0.05 in Rieke et al. (2023), and fout = 0 for 
42 sources at z > 8 in Hainline et al. (2023a). A range 
of potential redshifts z = 0.01 — 22 in Az = 0.01 in- 
crements were considered, and for selection, we adopt 
the use of the redshift corresponding to the minimum 
x? from the fit, za. For each nominal redshift, we use 
the Inoue et al. (2014) model for attenuation from the 
intergalactic medium (see also Madau 1995). We do not 
adopt any magnitude priors, we impose an error floor 
of 5% on the photometry, and allow for negative fluxes. 
When fitting the SED models to determine a photo- 
metric redshift, to maximize signal-to-noise ratios we 
use aperture-corrected r — 0.1" circular aperture fluxes 
on the native resolution JOF images without convolu- 
tion to a common PSF, multiplied by the photometric 
offsets discussed in Section 2.3. We have checked that 
we obtain comparably high-redshift solutions when us- 
ing common-PSF Kron aperture photometry with lower 
SNR, except where noted below. We note that for some 
objects, the best-fit SED model has Lyman-a line emis- 
sion. This feature arises as an artifact of the optimiza- 
tion process in EAZY that mixes templates with and 
without Lyman-a emission. We do not claim this line 
emission to be real. The equivalent width of Lyman-a is 
degenerate with the redshift of the break, which can con- 
tribute to a photometric redshift offset of Az & 0.2 — 0.4 
relative to a spectroscopic redshift. Local attenuation 
from the galaxy interstellar medium or circumgalactic 
medium can shift the photometric redshift by a similar 
amount (e.g., D'Eugenio et al. 2023; Heintz et al. 2023) 


3.2. Selection Criteria 


In the JOF, we apply the following criteria to iden- 
tify our high-redshift sample. These criteria have been 
adapted from Hainline et al. (2023a) but further tailored 
toa 12 S z < 20 selection. We note that these crite- 
ria both select objects previously discovered, notably by 
Hainline et al. (2023a), and identify new objects. We 
provide the provenance of each object when discussing 
our samples below. Our selection criteria are: 


1. The redshift at the EAZY fit y? minimum must be 
Za 2 11.5. 


2. Two of F277W, F356W, and F444W JWST NIR- 
cam filters must show > 5o detections. 


3. All the long-wavelength NIRCam fluxes (F250M, 
F277W, F300M, F335M, F356W, F410M, F444W) 
must exceed 1.50 significance. 


4. The redshift posterior distribution must have an 
integral probability of P(z, > 11) > 0.68, where 
we take P(z) x exp(—x?/2). 


5. The goodness of fit difference between the best 
high redshift (z > 11) and low redshift (z < 7) 
solutions must satisfy Ay? > 4, and for the best 
fit we require x? < 100 summed over all 19 filters. 


6. The flux in F090W and F115W each must be be- 
low 2.50 significance, as we expect no robust de- 
tection of flux blueward of the Lyman-q break. 


7. lo avoid objects redder than the typically blue 
high-redshift objects (e.g., Topping et al. 2023), 
we require that sources cannot have both (F277W- 
F356W)> 0.125 and (F356W-F444W)> 0.25. 


8. Each object must have F150W, F162M, F182M, 
F210M, and F277W coverage. 'This criterion limits 
our survey area to the F162M JOF footprint. 


9. The NIRCam short and long wavelength local ex- 
posure time must be within a factor of four, which 
avoids edge effects from the mosaic pattern. 


10. To avoid variable sources, the flux in the NIR- 
Cam Medium bands acquired in the second year 
of JWST operations must not exceed the broad 
band NIRCam fluxes acquired in the first year by 
more than 1c in all bands simultaneously. In prac- 
tice, we treat overlapping medium and broad band 
filters as random samples of the same flux density, 
and then flag when the difference between such 


pairs of flux estimates exceeds the quadrature sum 
of each pair's errors when taken in different epochs. 


11. We require that the source not be covered by an- 
other galaxy as determined from the segmentation 
map, which lowers the available area by 22%. The 
final available area after accounting for foreground 
sources is approximately A' — 7.06 square arcmin. 


We note that without the data quality (criteria 9- 
10), minimum long-wavelength SNR. threshold (cri- 
teria 2-3), or color criteria (criteria 6-7), fifteen 
objects would be selected. However, of these 
sources, one (JADES+53.05101-27.89787) sits in an 
oversubtracted area of a distant star diffraction 
spike and three more are covered by a stray light 
“wisp” feature in F162M (JADES+53.08317-27.86572, 


JADES+53.07681-27.86286, and JADES+53.04964- 
27.88605). For a discussion of wisp features in JWST, 
please see Rigby et al. (e.g., 2023). 


Of the remaining eleven, one fails the minimum SNR 
threshold (for JADES4-53.07385-27.86072, all filters 
redward of F335M have f, « 2nJy) and one fails the 
color criteria (JADES+53.08468-27.86666 is red). In 
total nine objects pass the selection; these comprise our 
Main Sample (see Table 2). While we will comment on 
the additional two interesting objects that nearly satisfy 
our selection, we do not consider them in our fiducial 
luminosity function analyses. We call this collection of 
two objects the “Auxiliary Sample" at z > 11.5. There 
are also five sources in Hainline et al. (2023a) sample in 
the vicinity of the JOF with previously reported pho- 
tometric redshifts z > 11 that are not in our sample. 
Of these, with the additional JOF data we find four 
objects to have revised photometric redshifts z « 7 or 
fail other selection criteria (JADES+53.02700-27.89808, 
JADES+53.03696-27.89422, JADES+53.07901- 
27.87154, JADES+53.10469-27.86187). The fifth falls in 
the F162M gap (JADES+53.07076-27.86544; NIRCam 
ID 176151) and therefore does not reside in our effective 
area. 

We note that the F250M SNR criterion fixes the upper 
redshift limit of our selection. If we remove this criterion 
and the z > 11.5 limit, we find one additional z ~ 11.4 
candidate (JADES+53.10131-27.85696) detected in all 
filters F150W and redder with f, ~ 2nJy, excepting 
F250M which is about 1.50 low. In other words, we 
would find no z = 20 candidates by removing the weak 
F250M detection criterion. 

The luminosity function analysis discussed below in 
850 enables the accounting of potential contributions 
to the inferred galaxy abundance from galaxies with 
maximum-likelihood photometric redshifts below the 
putative redshift of interest. We identified galaxies with 
maximum likelihood redshifts z > 8 and P(z > 12) > 
0.01 that otherwise satisfy the above selection criteria. 
There are three such galaxies, which fall in the pho- 
tometric redshift z ~ 10.5 — 11.2 range, which will be 
referred to as the “Contributing Sample" at z « 11.5. 
Two of these (NIRCam IDs 76035 and 172510) were pre- 
viously found in Hainline et al. (2023a). A third galaxy, 
NIRCam ID 64312 with photometric redshift z & 10.6 
and P(z > 12) ~ 0.05 from photometry, was subse- 
quently confirmed to lie at slightly lower redshift with 
P(z < 12) < 0.01 and was not considered further. 


Table 2. High-redshift candidates in the JADES Origins Field 


Name NIRCam ID RA Dec T Muv Rıj2 [mas]? ^ P(z < 7)° 
Main Sample z > 11.5 

JADES+53.09731-27.84714 74977 53.09731  -27.84714 11.53*027 —17.66 + 0.14 De 3.72% 10° 
JADES+53.02618-27.88716 16699 53.02618 -27.88716 11.56103} —17.94 + 0.15 351 9.12 x 1074 
JADES+53.04017-27.87603 33309 53.04017 -27.87603 12.10107 —17.73 +0.10 JA 4.02 x 10-5 
JADES+53.03547-27.90037 160071 53.03547 -27.90037 12.38*'01' —18.16+0.11 d. 7.87 x 1074 
JADES+53.06475-27.89024 13731 53.06475 -27.89024 12.93*0 95 —18.78 + 0.04 Cas Berg 
JADES+53.02868-27.89301 11457 53.02868 -27.89301 13.521026  —18.55 +0.11 1944 7.75 30-5 
JADES+53.07557-27.87268 376946 53.07557 -27.87268 1438753; —18.28 + 0.22 6*9 7.63 x 107? 
JADES+53.08294-27.855634 183348 53.082904 -27.85563  14.39*0 25 —21.00 + 0.05 HQE 0 

JADES+53.10762-27.86013 55733 53.10762 -27.86013 14.637029 —18.54 + 0.13 agt 2.26 x 107? 

Contributing Sample z < 11.5 
JADES+53.03139-27.87219 172510 53.03139 -27.87219 10.76Ł09$6 —17.85 + 0.10 827. 6.49 x 1075 
JADES+53.09292-27.84607 76035 53.09292 -27.84607 11.051039 —17.83 + 0.15 ES 4.06 x 1074 
JADES+53.06857-27.85093 70836 53.06857  -27.85093 11.17*029 —18.02 + 0.10 ME 2.38 x 1073 
Auxiliary Sample z > 11.5 

JADES+53.08468-27.86666° 44962 53.08468 -27.86666 12.9552 —18.16 + 0.10 561? 1.09 x 107? 
JADES+53.07385-27.86072/ 54586 53.07385 -27.86072 13.0670 25 —17.08 + 0.12 10° 6.16 x 10? 


NoTE—^ The half-light size refers to the intrinsic, PSF-deconvolved size of each source, in milliarcseconds. P Best fit photometric 
redshift with 16- and 84-percentile uncertainties from the inferred photometric redshift distribution. © The posterior probability 
density for the photometric redshift of the candidate to lie at redshift z « 7, given the SED fitting method described in 83.1. 4 
Spectroscopically confirmed at z = 14.32 by Carniani et al. (submitted). © Fails red color limit. f Fails minimum SNR criterion. 


Table 3. Aperture-corrected HST/ACS photometry? in r — 0.1" circular apertures. 


Name F435W [nJy] F606W [nJy] F775W |nJy] F814W [nJy]  F850LP [nJy] 
Main Sample z > 11.5 
JADES+53.09731-27.84714  —1.06 + 2.49  —0.58 + 1.35 2.20 + 3.52 —3.58 + 1.40 —11.28 + 4.63 
JADES+53.02618-27.88716 | —2.47 +3.10 —0.79 +3.71 —5.52 +5.67 2.57 + 2.33 7.01 + 7.62 
JADES+53.04017-27.87603 . 0.15 1.70  —1.38 +3.37 8.59 + 5.73 2.04 + 1.98 —2.17 + 7.90 
JADES4-53.03547-27.90037 | —3.50 + 1.73 —2.86 +2.40 —2.08 +3.28 1.55 + 2.24 8.51 + 8.26 
JADES+53.06475-27.89024 1.60 + 2.97 0.47 + 1.79 —3.36+2.67 4.45 + 2.30 0.20 + 5.04 
JADES4-53.02868-27.89301 —4.05 + 3.02 —1.82 +3.33 4.03 + 3.64 2.84 + 2.22 3.77 + 8.78 
JADES+53.07557-27.87268 | —0.99 +1.88 1.83 + 1.73 0.99 + 3.72 0.73 + 1.80 3.17 + 4.94 
JADES4-53.08294-27.85563 —2.80 +3.47 0.54 + 1.36 3.87 + 3.94 3.67 + 1.47 1.66 + 4.45 
JADES4-53.10762-27.86013 —3.64 + 2.62 —0.29 + 1.56 3.47 +3.38 —1.35 + 1.89  —0.35 + 4.72 
Contributing Sample z < 11.5 
JADES+53.03139-27.87219 —1.30 +1.69 3.82+485  —1.02 +5.67  —1.37 + 1.46 5.42 + 8.00 
JADES4-53.09292-27.84607 | —0.95 +2.41 0.69 + 1.27 0.76 + 3.56 | —2.21 1.41  —2.95 + 4.35 
JADES+53.06857-27.85093 4.54 + 2.64 2.21 + 1.33 —7.54+2.83 3.70 + 1.40 —1.87 + 3.75 
Auxiliary Sample z > 11.5 

JADES4-53.08468-27.86666 | —0.01 + 2.19 —1.82+1.37 3.5143.72  —0.03 +1.58  —6.87 + 4.91 
JADES+53.07385-27.86072 7.83 + 2.55  —0.31+1.38 0.98 +3.71 —0.17 + 1.34 0.43 + 4.35 


NoTE—* These photometric measurements were made on the native-resolution Hubble Legacy Field images 
(Illingworth et al. 2016b). 
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Figure 1. F444W/F200W/F090W false color red/green/blue image of the JADES Origin Field (background image; 27.5 
arcmin’), the JOE F162M footprint (jade outline) and F356W+F410M+F444W /F200W+F210M/F090W+F115W false color 
red/green/blue thumbnail images (each 0.86 arcsec”) for z > 12 high-redshift galaxy candidates. The RGB images of the galaxy 
candidates typically appear to have a green hue in this color space, as they are all detected in the filters used for both the 
green and red channels, but not the blue channel. Each inset thumbnail lists the best-fit EAZY photometric redshift and the 
JADES NIRCam ID, and we indicate the shared angular scale of the thumbnails with a scale bar showing 0.2". Table 2 lists 
the designations of the objects based on [RA, Dec]. NIRCam ID 183348 was spectroscopically confirmed as JADES-GS-z14-0 
by Carniani et al. (submitted) at z — 14.32. 


Table 4. Aperture-corrected short-wavelength JWST/NIRCam photometry in r = 0.1" circular apertures. 


Name FO9OW [nJy] F115W [nJy] F150W [nJy] F162M [nJy] F182M [nJy] F200W [nJy] F210M [nJy] 
Main Sample z > 11.5 
JADES+53.09731-27.84714 0.67 £0.53 —0.01 + 0.46 2.13 + 0.46 4.14 + 0.56 2.93 + 0.36 3.49 + 0.52 2.02 + 0.47 
JADES+53.02618-27.88716 0.62 + 0.61 —0.49 + 0.50 1.80 + 0.47 4.10 + 0.63 3.85 + 0.45 3.83 + 0.49 2.83 + 0.56 
JADES4-53.04017-27.87603  —0.20 + 0.59 —0.32 + 0.49 0.56 + 0.45 2.81 + 0.58 3.14 + 0.42 3.71 £0.45 3.43 + 0.49 
JADES+53.03547-27.90037 0.65 + 0.52 —0.69 + 0.42 0.65 + 0.41 2.06 + 0.53 3.06 + 0.38 2.49 + 0.41 2.82 + 0.48 
JADES4-53.06475-27.89024 | —0.08 + 0.50 0.29 + 0.40 0.10 + 0.37 1.16 + 0.46 7.80 + 0.36 6.24 + 0.42 7.48 + 0.41 
JADES4-53.02868-27.89301 —0.82 + 0.66 0.77 + 0.54 0.54 + 0.49 1.46 + 0.64 4.97 + 0.46 5.90 + 0.51 7.46 + 0.58 
JADES+53.07557-27.87268 0.21 + 0.52 0.21 + 0.43 —0.82 + 0.42 0.17 + 0.52 0.75 + 0.38 2.18 + 0.41 0.60 + 0.46 
JADES+53.08294-27.85563 —1.12 + 0.68 0.73 + 0.66 1.10 + 0.55 — 9.73 + 0.90 20.78 +0.58 29.66 + 1.14 
JADES+53.10762-27.86013 0.36 + 0.56 0.60 + 0.47 0.16 + 0.44 0.74 + 0.61 1.87 + 0.38 3.37 + 0.45 3.72 + 0.53 
Contributing Sample z < 11.5 
JADES+53.03139-27.87219 —0.41 + 0.67 0.03 + 0.56 3.07 + 0.52 4.70 + 0.65 4.02 + 0.47 3.70 + 0.56 2.95 + 0.57 
JADES+53.09292-27.84607 0.16+0.51 0.26 + 0.44 2.36 + 0.50 4.03 + 0.63 3.43 + 0.37 3.43 + 0.54 3.29 + 0.57 
JADES+53.06857-27.85093 0.49 + 0.48 0.16 + 0.39 2.36 + 0.40 3.95 + 0.49 3.70 + 0.31 3.06 + 0.39 3.86 + 0.43 
Auxiliary Sample z > 11.5 
JADES+53.08468-27.86666 0.05 + 0.48 —0.47 + 0.43 0.38 + 0.42 0.60 + 0.46 3.02 + 0.39 2.78 + 0.42 3.82 + 0.44 
JADES4-53.07385-27.86072  —0.52 + 0.46 0.26 + 0.40 0.59 + 0.38 —0.01 + 0.48 2.03 + 0.37 1.55 + 0.39 1.81 + 0.44 
candidates in our Main Sample, and the two auxiliary 
objects that have photometric redshifts z > 11.5 but 
that fail data quality or redness cuts. 
Next, in order of increasing photometric redshift, we 
introduce each galaxy candidate with some summary 
3.3. Sample discussion and a figure of the SED fits to the 0.1”-radius 


Given the selection criteria presented in 83.2, our en- 
tire 11.5 < z < 15 sample consists of eight galaxy candi- 
dates (our ^Main Sample"). Table 2 lists their designa- 
tions based on |RA, Dec], the internal JADES NIRCam 
ID, and the best-fit redshift za. Five of these objects (ID 
16699, 33309, 13731, 11457, and 55733) were previously 
identified in Hainline et al. (2023a); the other three are 
new here. We also record galaxy sizes measured from 
our Forcepho modeling in Table 2. For each object, we 
provide r — 0.1" circular aperture photometry for HST 
ACS bands in Table 3, and JWST 0.1”-radius circular 
aperture photometry for the NIRCam short-wavelength 
and long-wavelength filters appear in Tables 4 and 5, 
respectively. Note the fluxes we report in Tables 3-5 are 
measured on the native resolution images and are not 
convolved to a common PSF. 

Figure 1 shows an  F444W/F200W/F090W 
red/green/blue false color mosaic of the JOF re- 
gion. Of the 27.5 arcmin? area shown, about 9.05 
arcmin? has acceptable F162M coverage. The inset 
thumbnail images for each galaxy candidate show 0.86 
arcsec? regions with red/green/blue colors provided by 


F356W+F410M+F444W /F200W 4-F210M/F090W --F115W, 


along with the best fit redshift z, and the JADES NIR- 
Cam ID referenced in Table 2. We plot the eight galaxy 


circular aperture photometry. We show the photomet- 
ric redshift posterior distribution and best-fit redshift 
for each and the redshift posterior distribution limited 
to z « 7, as well as the best-fit SED and most likely 
low-redshift SEDs. We also show the JWST filter trans- 
mission curves and the fourteen JWST filter cutouts for 
each galaxy. 


3.3.1. JADES+53.09731-27.84714; NIR Cam ID 74977 


Figure 2 shows the best-fit SED for object 
JADES4-53.09731-27.84714 (NIRCam ID 74977). The 
object is remarkably faint with map œ~ 30.5 redward of 
Lyman-o, and has a best-fit redshift of z, = 11.5. The 
best low-redshift solution has Złow = 2.7, but exceeds 
the observed F115W constraint. 


3.3.2. JADES+53.02618-27.88716; NIR Cam ID 16699 


Figure 3 shows the best-fit SED for object 
JADES+53.02618-27.88716 (NIRCam ID 16699, Hain- 
line et al. 2023a). The best-fit redshift is z, = 11.6 
for this faint source, which has mag œ 30.2 — 30.5 in 
the NIRCam long-wavelength channels. The best low- 
redshift solution has zou = 2.6. We note that using the 
BAGPIPES SED-fitting code (Carnall et al. 2018) to 
constrain the photometric redshift of this galaxy candi- 


10 


Table 5. Aperture-corrected long-wavelength JWST/NIRCam photometry in r = 0.1" circular apertures. 


Name F250M |nJy] F277W [nJy] F300M [nJy] F335M [nJy] F356W [nJy] F410M [nJy] F444W [nJy] 
Main Sample z > 11.5 
JADES+53.09731-27.84714 . 1.68 + 0.73 2.13 + 0.33 2.49 d- 0.52 2.25 + 0.45 2.04 + 0.35 2.96 + 0.60 3.07 + 0.49 
JADES+53.02618-27.88716 1.74 + 0.66 2.15 + 0.35 2.32 + 0.46 2.57 + 0.36 2.17 + 0.38 1.76 + 0.63 2.50 + 0.51 
JADES+53.04017-27.87603 2.81 + 0.75 2.63 + 0.34 2.70 + 0.47 2.80 + 0.37 2.52 + 0.37 1.71 + 0.63 1.20 + 0.51 
JADES+53.03547-27.90037 . 3.16 + 0.65 3.56 + 0.34 3.58 + 0.44 2.73 + 0.38 2.32 + 0.39 2.79 + 0.66 2.47 + 0.52 
JADES+53.06475-27.89024 7.22 + 0.65 5.34 + 0.30 5.20 + 0.46 6.21: 0.37 5.47 + 0.34 5.01 + 0.58 4.30 + 0.46 
JADES+53.02868-27.89301 5.63 + 0.69 4.71 + 0.36 4.55 + 0.46 DU = Vor 5.50 + 0.38 5.79 + 0.65 5.24 + 0.51 
JADES+53.07557-27.87268 . 1.80 + 0.60 2.22 + 0.33 2.77 + 0.45 2.27 + 0.39 2.41 + 0.35 2.18 + 0.57 2.20 + 0.46 


JADES4-53.08294-27.85563 32.02 +0.75 32.83 +0.44 31.30 +0.55 27.43 +0.48 28.21 +0.44 28.56 +0.71 29.58+0.55 


JADES+53.10762-27.86013 3.90 Æ 0.70 3.75 + 0.38 4.25 + 0.50 4.36 + 0.45 3.76 + 0.42 2.68 + 0.65 4.07 + 0.49 
Contributing Sample z < 11.5 
JADES+53.03139-27.87219 . 2.34 + 0.69 2.(2-:5 0,93 3.51 + 0.47 2.90 + 0.38 2-20 0.50 2.21 + 0.62 3.33 + 0.51 
JADES--53.09292-27.84607 3.96+0.71 2.41 + 0.34 3.26 + 0.52 2.51 + 0.46 2.43 + 0.37 2.90 + 0.60 1.69 + 0.49 
JADES+53.06857-27.85093 . 4.49 + 0.68 3.63 + 0.37 3.08 + 0.49 2.87 + 0.43 2.98 + 0.37 3.21 + 0.59 2.25 + 0.48 
Auxiliary Sample z > 11.5 
JADES+53.08468-27.86666 — 4.31 + 0.73 4.21 + 0.36 5.08 + 0.50 5.14 + 0.45 4.78 + 0.39 5.30 + 0.62 6.14+0.51 
JADES+53.07385-27.86072 1.62 + 0.69 2.12 + 0.39 2.60 + 0.51 2.28 + 0.44 1.84 + 0.39 1.41 + 0.64 1.74 0.51 


date provides a slightly lower redshift of z ~ 11.3, with- 
out Lyman-o emission and with still non-zero P(z » 
11). 


3.9.9. JADES+53.04017-27.87608; NIR Cam ID 33309 


Figure 4 shows the best-fit SED for object 
JADES4-53.04017-27.87603 (NIRCam ID 33309, Hain- 
line et al. 2023a). The best-fit SED model has za = 12.1, 
while the best low-redshift solution has zj5,, = 3.2. The 
source is also remarkably faint, with mag œ% 30.2 in the 
NIRCam long-wavelength filters. 


3.9.4. JADES+53.03547-27.90037; NIRCam ID 160071 


Figure 5 shows the best-fit SED for object 
JADES+53.03547-27.90037 (NIRcam ID 160071). The 
flux densities of this object are f, ~ 3.5 nJy (map & 
30). When fitting an SED model, the observed photom- 
etry, including the strong break in F150W, constrain 
the redshift to be z, — 12.4. The best solution at low 
redshift has Zow = 3.4. 


3.3.0. JADES+53.06475-27.89024; NIRCam ID 13731 


Figure 6 shows the best-fit SED for object 
JADES+53.06475-27.89024 (NIRCam ID 13731, Hain- 
line et al. 2023a). The long-wavelength JWST/NIRCam 
photometry shows m Ap ~ 29.5, and constrains the pos- 
terior photometric redshift distribution to be peaked 
strongly near z, = 12.9. The best-fit low-redshift so- 
lution at Zow = 3.5 would exceed the F090W, F115W, 
F150W, and F162M photometry by several standard de- 
viations. 


3.3.6. JADES+53.02868-27.89301; NIRCam ID 11457 


Figure 7 shows the best-fit SED for object 
JADES4-53.02868-27.89301 (NIRCam ID 11457, Hain- 
line et al. 2023a). The object has NIRCam flux densi- 
ties redward of the break of f, ~ 4 — 6 nJy (map & 
29.5 — 29.9), that constrain SED models to yield a best- 
fit redshift z, — 13.5. The best-fit low-redshift solu- 
tion SED at zou = 3.6 exceeds the F090W and F150W 
constraints by several standard deviations. The second 
peak in the high-redshift p(z) at z & 12.7 is driven by 
the marginal (2.30) detection in F162M, which if real 
would prefer a slightly lower redshift than the mode but 
still within our selection criteria. However, we caution 
that the F162M detection, the F182M-F210M color, and 
the rising SED shape longward of 3.5um could indicate 
a potential low-redshift contaminant not well-modeled 
by our SED template set. We therefore proceed with 
caution while including this candidate in our sample. 


3.9.(. JADES+53.07557-27.87268; NIR Cam ID 876946 


Figure 8 shows the best-fit SED for object 
JADES+53.07557-27.87268 (NIRcam ID 376946). This 
faint (map = 30.5) object at redshift za = 14.4 
is slightly redder than most of the other candidates. 
JADES+53.07557-27.87268 displays an unusual SED in 
that either the F182M and F210M fluxes must be biased 
low by several sigma to be consistent with the F200W 
flux, and the high-redshift solution does not match well 
the observed F182M, F200W, and F210M data. The 
best solution at low redshift has Zlow = 3.8 with nearly 
896 of the EAZY probability, although it overpredicts 
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Figure 2. SED model, photometric redshift posterior distributions, and JWST NIRCam image thumbnails for galaxy candidate 
JADES+53.09731-27.84714 (NIRcam ID 74977). The upper left panel shows the aperture-corrected r = 0.1" flux density f, in 
nJy of the NIRCam (purple points with lo uncertainties) and HST/ACS (red points with 1e uncertainties) photometry for the 
object, with median photometric offset corrections applied. The best-fit SED is shown in blue, while the best fit low-redshift 
solution is shown in gray. The synthetic model photometry for both models are shown as open squares, and the JWST NIRCam 
filter transmission curves are shown as colored regions. The upper right panel shows the posterior distribution of photometric 
redshifts for the object (blue), the best-fit redshift (vertical dashed line), the photo-z posterior if only redshifts z « 7 are 
considered (light gray), and the best-fit redshifts provided as an annotation, as is the posterior probability density at redshifts 
below z ~ 7. The bottom panel shows inverted grayscale thumbnails of the fourteen NIRCam filters in a 0.93 x 0.93 arcsec? 
region around each object, the stretch applied to each filter scaled with the mean value in the thumbnail. The signal-to-noise 
ratio of the aperture-corrected r — 0.1" circular aperture photometry for each band is noted in the corresponding thumbnail. 


The JADES NIRCam ID is also provided on the left side of the image. 


the observed F150W flux. We note that when using 
the BAGPIPES SED-fitting code (Carnall et al. 2018) 
with a broad log-uniform prior (M, € [10?,10!?] M5) 
on stellar mass to constrain the photometric redshift of 
this galaxy candidate, we find a yet larger low-redshift 
probability density than with EAZY. The best fit red- 
shift is still z > 14 and most of its photometric redshift 
posterior probability is at very high-redshift. We also 
note that this object has the largest increase in the low- 
redshift probability density when using common-PSF 
Kron aperture fluxes to fit a photometric redshift, but, 
given the loss in SNR for this exceedingly faint object, 
the photometric SED become much noisier. 


3.3.8. JADES+53.08294-27.8556; NIRCam ID 183348 


JADES4-53.08294-27.8556 (NIRcam ID 183348) with 
redshift z — 14.4 is the most remarkable object in our 
sample, with a best-fit SED shown in Figure 10. The 
object appears relatively bright (f, ~ 30nJy; r = 0.1” 


radius aperture) but shows strong break from F210M 
to F182M and no significant flux at shorter wave- 
lengths. Before the JOF ultradeep JWST/NIRCam 
medium band data was acquired, based on JADES 
JWST/NIRCam broadband data Hainline et al. (2023a) 
first discussed this source with a photometric redshift of 
Zphot — 14.51. Owing to the observed brightness of the 
source and its close proximity to another lower-redshift 
source, 183348 was rejected from their main sample. 
Subsequently, Williams et al. (2023) determined a lower 
photometric redshift z,not = 3.38, and found the source 
was detected by JWST/MIRI at 7um from the SMILES 
program (PID 1207; PI Rieke). Given the addition of 
our ultradeep JOF JWST/NIRCam medium band data, 
we find the photometric redshift posterior of 183348 dis- 
tribution is sharply peaked at z ~ 14.4. This high red- 
shift peak is now much more strongly favored than low 
redshift solutions as the new JOF medium band mea- 
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Figure 4. Same as Figure 2, but for galaxy candidate JADES+53.04017-27.87603 (NIRCam ID 33309). 
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Figure 5. Same as Figure 2, but for galaxy candidate JADES+53.03547-27.90037 (NIRCam ID 160071). 
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Figure 6. Same as Figure 2, but for galaxy candidate JADES+53.06475-27.89024 (NIRCam ID 13731). 
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Figure 7. Same as Figure 2, but for galaxy candidate JADES+53.02868-27.89301 (NIRCam ID 11457). 
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Figure 8. Same as Figure 2, but for galaxy candidate JADES+53.07557-27.87268 (NIRCam ID 376946). 
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Figure 9. Same as Figure 2, but for galaxy candidate JADES+53.10762-27.86013 (NIRCam ID 55733). 
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Figure 10. Same as Figure 2, but for galaxy JADES+53.08294-27.85563 (NIRCam ID 183348). We note this object has been 
discussed previously in Hainline et al. (2023a) and Williams et al. (2023), and spectroscopically confirmed by Carniari et al. 
(submitted). The F162M data for this object has been omitted because of data quality issues. 
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surements better constrain the shape and depth of the 
break at ~ 1.8um while placing limits on strong emission 
lines redward of the break. While low-redshift solutions 
have low probability, the low-redshift photometric red- 
shift posterior distribution is very sharply peaked near 
Zlow = 3.4 and requires a very red object with strong 
emission lines in F200W and F277W. A principal con- 
cern regarding 183348 is the close proximity of a neigh- 
boring galaxy (NIRCam ID 183349) that has a best-fit 
photometric redshift of za & 3.4. This alignment obvi- 
ously supported the previous suspicion that 183348 was 
also at the lower redshift. However, our analysis of the 
initial JOF NIRCam medium-band photometry as well 
as JWST/MIRI photometry (Helton et al., submitted) 
further supported the higher redshift, and on that ba- 
sis, the galaxy was selected for spectroscopic followup. 
Carniani et al. (submitted) present a spectroscopic red- 
shift confirmation of z = 14.32, and we refer the reader 
to that work for a detailed analysis of the properties of 
this intriguing galaxy. Here, we do compare the prop- 
erties inferred for this galaxy along with other objects 
in the Main Sample measured in the same manner. We 
note that the photometric and spectroscopic redshift dis- 
tributions are very similar, and our choice to adopt its 
photometric redshift distribution during the luminosity 
function inference has little impact on our results. We 
also note that gravitational lensing by the neighbor is 
considered by Carniani et al. (submitted), but find the 
magnification to be small. 


3.9.9. JADES-+53.10762-27.86013; NIRCam ID 55733 


Figure 9 shows the best-fit SED for object 
JADES4-53.10762-27.86013 (NIRCam ID 55733, Hain- 
line et al. 2023a). The galaxy candidate has NIRCam 
long-wavelength fluxes of map ~ 29.9 and a best-fit 
redshift of za ~ 14.6. The best low-redshift solution 
has Zow = 3.9 with 2% of the EAZY probability, al- 
though the corresponding SED model would substan- 
tially exceed the observed F150W. We note that this 
object shows F162M flux at 1.10 significance, and con- 
firmation of this hint of a signal would negate a possible 
high-redshift solution. 


3.9.10. Auxiliary Objects 


We also provide SED fits for the Auxiliary can- 
didates JADES+53.07385-27.86072 (Figure 11), and 
JADES+53.08468-27.86666. (Figure 12). 

JADES+53.07385-27.86072 (NIRcam ID 54586) is ex- 
ceedingly faint and is relegated to our Auxiliary sam- 
ple by failing the minimum SNR criteria, with some 
long-wavelength NIRCam filters showing map > 30.5 
flux levels. The high-redshift posterior distribution for 


this object is correspondingly broader, with a peak at 
Za = 13.1. The best low redshift solution has Zow = 3.6. 

Finally, JADES4-53.08468-27.86666 (NIRCam ID 
44962, Hainline et al. 2023a) is in our Auxiliary sample 
owing to its red SED that increases from f, 3 nJy in 
F182M to f, ~ 6nJy in F444W. The redshift posterior 
distribution is double-valued, with a peak at z, = 12.9. 
The best low-redshift solution has Zlow = 3.5. 


4. COMPLETENESS SIMULATIONS 


The detection and selection of high-redshift galaxy 
candidates impose limitations that reduce the complete- 
ness of a sample. To convert the number of observed 
galaxies satisfying the selection criteria into a measure- 
ment of the galaxy number density, the completeness 
of the detection and selection process can be computed 
and incorporated. Below, in 84.1 we use simulations to 
characterize our detection completeness and in 84.2 we 
simulate our selection completeness. These calculations 
are used in $5 to include completeness corrections in the 
rest-UV luminosity function. 

We note that the requirement to compute the com- 
pleteness suggests that the detection and selection pro- 
cess should be algorithmic and automatable. We there- 
fore do not apply any cuts based on visual inspection 
or judgment beyond crafting the detection method de- 
scribed in 82.3 or the selection criteria presented in $83. 
This restriction allows us to simulate both the detection 
and selection completeness. 


4.1. Detection Completeness 


To compute the detection completeness of our photo- 
metric pipeline, we performed detailed source injection 
simulations using a wide range of input sources. First, 
we create a mock input galaxy catalog by drawing from 
randomized distributions of galaxy physical properties 
including redshift, star formation rate, stellar mass, size, 
Sérsic (1968) surface brightness profile index, position 
angle, and axis ratio. The objects are selected to have 
properties comparable to the z > 8 sources reported by 
Hainline et al. (2023a). We use the Prospector code 
(Johnson et al. 2021) to compute the object fluxes given 
their physical properties and redshift. With this mock 
catalog, we use the GalSim (Rowe et al. 2015) image sim- 
ulation software to create simulated Sérsic (1968) profile 
objects distributed across a grid on the sky. We compute 
the overlap of the JOF mosaics in each filter with this 
grid of objects, and then add the randomized objects as 
injected sources in the JOF images. The result is a large 
set of synthetic JOF mosaics with injected sources. We 
can then process the images identically to the real data 
and attempt to discover sources. 
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Figure 12. Same as Figure 2, but for Auxiliary galaxy candidate JADES+53.08468-27.86666 (NIRCam ID 44962). 
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With the injected images, we combine the long- 
wavelength NIRCam images as for the real data, creat- 
ing an ultradeep stack. Our pipeline detection algorithm 
is applied to the injected mosaic stack to create a new 
detection catalog with simulated sources. We can then 
characterize the completeness of our detection method 
as a function of the source properties. We repeat the 
simulations with ten separate realizations, such that a 
total of 115,000 injected sources with widely-ranging in- 
trinsic properties are used. 


Figure 13 shows the detection completeness as a func- 
tion of the two main factors affecting this completeness. 
The apparent brightness of the objects influence their 
signal-to-noise ratio in the stacked detection image. The 
size of the object affects the surface brightness, which 
in turn determines the per pixel SNR that governs the 
contrast an object of a given luminosity relative to the 
sky background. The detection algorithm reaches 9096 
completeness at around map ~ 30.2 for small objects 
(Ri/2 S; 0.1 arcsec). This completeness function can be 
integrated into an interpolator to allow for the object 
completeness as a function of apparent magnitude and 
size to be utilized in inferring the UV luminosity func- 
tion. We note that through this simulation for the JOF 
we find that only about 78% of the pixels are not im- 
pacted by foreground objects, which we account for in 
computing our effective survey volume. Given that the 
objects of interest are small, only several pixels across, 
and our detection method reaches fairly low significance 
(SNR~ 1.5) per pixel such that the segmentations reach 
low surface brightnesses, we find this number to be rep- 
resentative of the impact of foreground sources on our 
detection completeness. 


4.2. Selection Completeness 


To simulate the selection completeness, we can use 
the spectral energy distributions in our mock galaxy 
catalog and the photometric uncertainty measured for 
our galaxy sample to simulate the effects of photomet- 
ric noise on our selection and consequently the inferred 
UV luminosity function. We create a sample of two mil- 
lion mock galaxies with model SEDs, induce photomet- 
ric noise with a normal scatter in each HST and JWST 
filter of the magnitude of our measured sky background. 
Our measurement uncertainties are sky-dominated, so 
only include sky noise in our simulated fluxes. These 
two million noisy simulated SEDs are then provided to 
EAZY exactly in the same manner as our real catalog, and 
SED fitting is performed to each object. This enables us 
to estimate how the photometric noise can disrupt the 
mapping between true redshift and photometric redshift, 
and identify which redshift windows could provide non- 


Intrinsic Size log mur: [arcsec] 
Completeness 
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Figure 13. Detection completeness in our JOF analysis as 
a function of intrinsic half-light radius and F277W apparent 
magnitude. The detection method is complete for small ob- 
jects and bright magnitudes, and the differential complete- 
ness reaches about 90% at F277W% 30.2AB for small ob- 
jects. Shown is a two-dimensional normalized histogram of 
object size and flux indicating the fraction of sources with 
such properties detected by the pipeline. The method be- 
comes highly incomplete fainter than map ~ 31 or for half- 
light radii above about half an arcsec. Owing to pixels cov- 
ered by foreground sources, the maximum detection com- 
pleteness will be reduced to ~ 78% of that shown here. 


negligible contamination for our selection criteria. For 
reference, we note that in our simulations, the fraction 
of objects with F200W SNR> 5 that are photometric 
redshift outliers with (|za — Ztrue|/(1 + Ztrue)) > 0.1 is 
3.8%. 

Figure 14 shows the completeness of selection crite- 
ria as applied to our mock galaxy catalog, as a function 
of the true object redshift and absolute UV magnitude. 
The selection proves highly complete at Myy < —18 for 
redshifts z = 12. At magnitudes fainter than Myy > 
—17.5, the photometric noise prevents the strict elimi- 
nation of low-redshift solutions such that the objects fail 
the Ax? selection described in 3. At the high-redshift 
end, the selection declines at z ~ 20 when the Lyman-a 
break affects F250M and our SNR requirement in that 
filter becomes limiting. As with the detection complete- 
ness, an interpolator can be constructed from the selec- 
tion completeness and then used to correct the galaxy 
number counts for the lossy selection process. We note 
here that we define Myy as the rest-frame 1500À UV 
luminosity density, computed by fitting a power-law to 
rest-frame UV photometry and marginalizing over any 
covariance with the spectral slope (for more details, see 
86.1). 
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Figure 14. Completeness of our selection criteria as a func- 
tion of galaxy redshift and absolute magnitude. For bright 
objects, the selection criteria described in 83 produce a sub- 
stantially complete sample. For fainter objects, the Ax? cri- 
terion fails as the photometric noise prevents the SED fitting 
procedure from distinguishing robustly between high and low 
photometric redshifts. 


5. REST-FRAME UV LUMINOSITY FUNCTION 
AT Z> 12 


To compute the rest-frame UV luminosity function 
from our sample of galaxy candidates and our complete- 
ness calculations, we can construct multiple measures of 
the galaxy abundance. We wish to account for several 
confounding effects. 

First, galaxies with a range of intrinsic redshifts will 
contribute to the the observed number counts of galax- 
ies at a given photometric redshift. The degree of this 
contamination will depend on the abundance of galaxies 
at other proximate intrinsic redshifts whose photomet- 
ric redshifts overlap with the epoch of our measurement. 
We must therefore account for the evolving luminosity 
function and mixing between populations at different 
redshifts. 

Second, each individual galaxy has a posterior distri- 
bution for its photometric redshift. Rather than assign 
each galaxy to a specific redshift bin and absolute mag- 
nitude, we can allow for a posterior distribution on the 
photometric redshift to represent a track of inferred ab- 
solute magnitude and redshift. Each galaxy can make 
a fractional contribution to the UV luminosity function 
at redshifts where its posterior has support. 

Given these considerations, we want to allow for flex- 
ibility in our representation of the UV luminosity func- 
tion. We can either infer a parameterized luminosity 
function by computing the likelihood of observing each 


19 


galaxy, given the evolving distribution of galaxy counts 
with luminosity and redshift, fully without binning, or 
we could bin in magnitude and redshift but account for 
the photometric redshift posterior distributions of each 
object. In either case, with the known individual proper- 
ties of each object, we want to treat the completeness of 
our detection and selection methods at the per-object 
level rather than through binning. Below, we present 
both methods, where we expand on the methods used 
by Leja et al. (2020) to infer the evolving stellar mass 
function at low redshift but now applied to the UV lu- 
minosity function evolution at high redshifts. We have 
tested both methods using mock galaxy samples con- 
structed from specified luminosity functions and poste- 
rior photometric redshift distributions. 


5.1. Inferring Evolving Luminosity Function 
Parameters 


The probability of observing an object with a given 
true luminosity and redshift is given by the product 
of the redshift dependent luminosity function (L, z|0), 
the selection function S(L,z), and the differential co- 
moving volume element probed V(z). We can assume 
the luminosity function depends on some parameters 0. 
Unfortunately, we do not know the true luminosity and 
redshift of each galaxy i, but instead estimate it from 
photometric data D;, by using SED models to construct 
the likelihood function £(D;|L,z). The likelihood of 
observing a galaxy with D; must then be marginalized 
over the unknown true parameters 


E (Dj | 0) x f dL | dz £(Di|L,z)A(L,z|0) (1) 
ACL, z|8) = (L, z|0) S(L, z) V(z) (2) 
$(Lz8) = é(z)(D/L.(z2))* De ro (3) 


Here A(L,z) is the differential number of objects ex- 
pected to be selected from the survey, as a function of 
the true L and z. We have parameterized the luminos- 
ity function as a single Schechter function. The red- 
shift evolution of the luminosity function can be treated 
with a dependence of the parameters on (z — 2,5) where 
Zref ls some reference redshift, e.g. the midpoint of the 
redshift range of interest. For our purposes, we will 
adopt either simple log-linear or log-exponential evolu- 
tion with redshift. To compute the likelihood of each 
object marginalized over the true object redshift and 
luminosity we numerically integrate the marginalization 
integrals using samples from the probability distribution 
provided by EAZY. 


CU: |0) ~ 5j; wi, ALi js Zij 


0)/2,,wij (4 


By drawing fair samples from the probability distribu- 
tions provided by EAZY, and noting that the effective 
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priors on z and L were uniform, each sample has equal 
weight w; j. With the ability to compute the likelihood 
of each object given the model, the likelihood for an en- 
semble of objects is then the product of the individual 
likelihoods. However, we must include the overall con- 
straint given by the number of observed objects. The to- 
tal expected number of selected objects is given by the 
integral of the product of the luminosity function and 
the effective volume, and the observational constraint is 
given by the Poisson likelihood of the actual number of 
observed objects! 


L(D|0)= e^" [T; £(Di | 6) (5) 
No = f dL f dz X(L, 2) (6) 


Here Na is the total number of observed objects. Note 
that for redshifts and luminosities for which our ob- 
servations are complete, the method accounts for the 
likelihood of non-detections given the chosen luminosity 
function parameter values. 


5.2. Estimating a Step-Wise Luminosity Function 


While the method in $5.1 does not bin in redshift or 
luminosity, the observed candidate galaxies could be as- 
signed to specific redshift and luminosity bins. If noth- 
ing else, binning allows for the measured galaxy abun- 
dance to be usefully plotted and compared with other 
measurements. The binned luminosity function summa- 
rizes the information retained by the unbinned param- 
eterized LF for which representing constraints on the 
galaxy abundance requires access to samples of the pos- 
terior distribution of LF parameters. 

Consider the photometric redshift posterior distribu- 
tion p;(z) of a candidate galaxy 7 with observed apparent 
magnitude m;. In the absence of photometric noise, the 
absolute magnitude of the object is M; = m; — DM(z), 
where DM (z) is the cosmological distance modulus in- 
cluding K-correction. Accounting for photometric noise, 
we will instead have some distribution of absolute mag- 
nitudes p(M;|m,;,z) for each object at a given photo- 
metric redshift. 'The distribution of inferred absolute 
magnitudes in some redshift bin zı to 22 is 


p(Mi|21, 22) = | dz | amp Mmi. pe) (7) 


The contribution of a galaxy to an absolute magnitude 
bin would then be 


M» 
N;(Mj, Ma, 21,22) = | p(M;|21,22)dM;. (8) 
Mı 


! This can be derived from the treatment of the luminosity func- 
tion as an inhomogeneous Poisson process; in the case that the 
effective rate À is constant this reduces to the typical Poisson 
likelihood. 


The total number density per magnitude n; in a mag- 
nitude bin Mı < M; < M» would then be 
2 4 NM, M2, 21, 22) 

n5 (M1, Ma, 21, 22) (Mo - M)V, (9) 
where V; is the average effective volume in the bin, al- 
lowing for the completeness to vary for each object i. In 
practice, evaluating these equations involves summing 
over samples from the photometric posterior distribu- 
tions of the galaxies while accounting for samples that 
lie outside the redshift bin to enforce the posterior nor- 
malization constraint | p(z)dz = 1. We note that when 
computing the samples in Myv and z, to compute Myv 
we use the 1500A rest-frame flux computed in the ap- 
propriate JWST filter given a putative redshift z. When 
computing Myv, we use the total fluxes computed from 
the Forcepho morphological decompositions. 

Procedurally, for each redshift bin we take all ordered 
Myv samples and separate them into bins whose edges 
are set to maintain a comparable number of samples per 
bin. We sum the number of samples in each bin and di- 
vide by the total number of samples across all galaxies, 
which provides the (non-integer) number of galaxies per 
bin. The average completeness in the bin is computed 
from the per-object selection and detection complete- 
ness based on the object properties and the fraction of 
pixels in the image not covered by foreground sources. 
We then divide the number of galaxies in each bin by 
the bin width, the completeness, and the volume to get 
the number density. The uncertainties for each bin are 
estimated from number count statistics. 

While we report our step-wise estimate, which ac- 
counts for photometric scatter between magnitude bins 
and variable completeness, we consider these measure- 
ments estimated checks on the inferred LF constraints 
described in 85.1 that do not bin in either redshift or 
magnitude and additionally account for potential con- 
tamination from proximate redshifts and the evolving 
shape of the luminosity function with redshift. We em- 
phasize here that our formal derived constraints on the 
luminosity function are provided through our inference 
procedure in the form of the computed posterior distri- 
butions of the parameters of our model evolving lumi- 
nosity functions. 


5.3. Luminosity Function Constraints 


Given the measured properties of our sample galaxies, 
their photometric redshift distributions p(z), and the 
method described in 85.1, we can compute marginalized 
constraints of an evolving UV luminosity function once 
we adopt a parameterized form. 


Table 6. Step-wise Luminosity Func- 


tion. 
Muv duu [107^ mag! Mpc ^ ?] 

11.5 « z « 13.5 

—-18.5*5 15 1.22 + 0.94 

“18.0037. 3.20 + 2.46 

+0.65 

zd ur 1.54 + 1.18 
15515 

—20.8* 215 0.971 + 0.357 

—18.4*0 10 2.56 + 2.46 

-181*113 0.783 + 0.754 


NOTE— The ranges listed for each Myv re- 
flect the widths of the magnitude bins, 
which are determined by the distribution of 
photometric redshift posterior samples for 
the objects contributing to each bin. 


Table 7. Luminosity Function Marginalized Param- 
eter Constraints 


Parameter Prior Constraint 
login ?9  U(-8,-2) | —6.339  —5222  —4.24 
M* u(—17,—24) -—24.95 —22.80 —20.71 
n° U(—3, 3) —0.29  —0.20 —0.13 
a U(—3,—1) 2160 —LT79 4 


NoTE— % The lower limit on the LF normalization is not 
well constrained, but the 95% upper limit is log,0 dn < 
—3.84. © The 95% upper limit on the characteristic mag- 
nitude is M* « —19.9. * We constrain the evolution 
parameter to be 7 < —0.08 at 95%. 


For the luminosity function, we adopt a redshift- 
dependent Schechter (1976) function 


buv(Muv, z) = 0.4log 10¢,(z)[10° 404" "Movie"? 


where the redshift-dependent normalization $,(z) can 
be further parameterized. Our fiducial choice for the 
normalization evolution is 


log10 $,(2) = logio Px,0 + (2 — 20). (11) 


We will refer to z as the reference redshift, which 
we will take fixed at z = 12 unless otherwise noted. 
The default parameters of the model then include d = 
[M*, a, 6,0, |, or the characteristic magnitude M*, the 
faint-end slope a, the normalization at the reference red- 
shift o, o, and the log-linear rate of change with redshift 
n. In practice, we fit in maggies | = —0.4Myy and 
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then convert to absolute magnitudes after inference. We 
adopt log-uniform priors for o, o and 7, a uniform prior 
in magnitude, and a uniform prior in a. The priors are 
reported in Table 7, along with our inferred constraints 
on the parameters. We emphasize again that informa- 
tion from all redshifts where the selection function has 
non-negligible support is included by our inference pro- 
cedure, which accounts both for regions of redshift and 
magnitude space with detections and those absent sam- 
ples that could have been detected if present. ‘The ef- 
fective redshift range where our model is informative for 
the luminosity function is mostly set by the selection 
completeness (Figure 14), or roughly z ^ 11 — 20. We 
present the full posterior distributions on the parameters 
in Figure 15. We here emphasize that the clear covari- 
ance between o, and M* mostly acts to keep the lumi- 
nosity density pyy x L*o, roughly constant at a given 
redshift. This feature is reflected in our constraints on 
puv shown in Figure 17. 

Since we constrain the abundance of galaxies at all 
selected and detectable redshifts and magnitudes simul- 
taneously, evaluating the luminosity function at any one 
redshift requires computing the marginal distribution of 
the luminosity function equation 10 over the posterior 
distribution of parameters for a given redshift and range 
of absolute magnitudes. At each z and Myy, equation 
10 is evaluated for all posterior samples, and the cu- 
mulative distribution of pyy weighted by the sample 
weights ww; constructed. Figure 16 shows the marginal 
constraint on the UV luminosity function at redshift 
z = 12, with the 16-84% of rv shown as a shaded 
region and the median ópyy shown as a white line. We 
also show the median inferred @yy at z = 14 as a light 
gray line. Note that none of these yy percentiles are 
guaranteed to follow equation 10 individually, but we do 
report the marginalized constraints on the luminosity 
function parameters in Table 7. We also show our step- 
wise luminosity function estimates computed in redshift 
bins of 11.5 < z < 13.5 and 13.5 < z < 15. These step- 
wise luminosity function measures are reported in Table 
6. 

In Figure 16, we also show z ~ 12—14 luminosity func- 
tion determinations reported in the literature. ‘These 
measurements include the z ~ 12 data from Harikane 
et al. (2023b), Harikane et al. (2023a), Pérez-González 
et al. (2023a) and Willott et al. (2023), the Adams et al. 
(2023b) constraints at z ~ 12.5, z ~ 13 measurements 
from Donnan et al. (2023a) and McLeod et al. (2023), 
and the z ~ 14 determinations from Finkelstein et al. 
(2023a). 'The median luminosity function constraints 
inferred from our sample and our forward modeling ap- 
proach agree with the available observations to within 
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Figure 15. Posterior distributions of the evolving luminosity function parameters. Shown are the posterior distributions for 
the luminosity function normalization log,, d, [Mpc? mag !], the normalization evolution parameter 7, the characteristic 
magnitude M” in absolute magnitude, and the faint-end slope a. Contours represent the 68% and 90% enclosed probabilities 
for each parameter. The marginalized posterior distributions for each parameter are shown at the top of each column, along 
with the 16%, 50% and 84% marginal constraints (see also Table 7). The lower limits on $, and M, are not well constrained, 
but we constrain at 95% probability that log) 9. < —3.84 and M, < —19.9. We note that 7 < 0 with > 95% probability, 
indicating that we infer a declining luminosity density at z > 12. 


about lo, excepting the z ~ 14 constraints from Finkel- 
stein et al. (2023a) that lie above our inference. We note 
here that the z ~ 11 luminosity function constraints 
from Donnan et al. (2023a), McLeod et al. (2023), and 
Finkelstein et al. (2023a) lie above our 84% inference of 
the z = 12 luminosity function, and that our selection 
function (Figure 14) by design removes z ~ 11 galaxies 
from our sample. We also emphasize that our results 
are completely independent of the other data shown in 
Figure 16. 


9.9.1. Luminosity Density Evolution 


Given the evolving luminosity function parameters in- 
ferred given the sample properties, the UV luminos- 
ity density evolution pyy(z) can be computed. Fig- 
ure 17 presents the marginalized constraints on the UV 
luminosity density evolution. Shown are 16-84% (jade 
shaded region) and median pyy (white line) integrated 
to Muy < —17, along with measured (left panel) or 
extrapolated (right panel) constraints to Myy < —17 
from the literature. Our measurements have sensitivity 
to objects at redshifts 11 < z < 20, and we indicate the 
luminosity density evolution inferred for the model rep- 
resented by equations 10 and 11. As the Figure shows, 
we infer that the UV luminosity density declines at high- 
redshift at a rate of n = dlog¢,/dz ~ —0.2 per unit 
redshift. Between z = 12 and z = 14, we therefore 
infer that the luminosity density declines by a factor of 
10-9.2(14-12) ~ 2.5. Within our statistical uncertainties, 
this inference agrees with almost all the literature deter- 
minations including Ishigaki et al. (2018), Bouwens et al. 
(2022), McLeod et al. (2023), Donnan et al. (2023b), 
Harikane et al. (2023b), Harikane et al. (2023a), Adams 
et al. (2023b), Pérez-González et al. (2023a), Leung 
et al. (2023), and Willott et al. (2023). The constraints 
at z ~ 11 from Finkelstein et al. (2023a) agree with our 
results, but their z ~ 14 point lies above our constraints 
albeit with large uncertainties. If we extrapolate the 
UV luminosity evolution inferred by our model, we find 
good agreement with the literature measurements back 
to z ~ 8 (e.g., Ishigaki et al. 2018; Bouwens et al. 2022; 
Pérez-González et al. 2023a; Willott et al. 2023; Adams 
et al. 2023b). Also shown in Figure 17 is the correspond- 
ing evolution in the cosmic star formation rate density 


psrR, using the approximate conversion from puy of 
kuy = 1.15 x 10 75 Mo yr ! erg ! s Hz from Madau 
& Dickinson (2014). For comparison, we also show the 
Madau & Dickinson (2014) model for the evolving cos- 
mic star formation rate density. 


5.4. Caveats 


Of course, with only nine objects at these extreme dis- 
tances and depths, there are important caveats to con- 
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sider about the LF measurement. First, most of our ob- 
jects are photometric candidates, and despite the closer 
spacing of the medium bands and our care in selection, 
we consider it possible that some might be lower red- 
shift interlopers. A Lyman-o break at z = 14 falls at 
the same wavelength as a Balmer break around z & 4. 
We stress that false positives would likely have a redshift 
distribution that falls less slowly than the true Lyman-a 
break population, so a population of false positives will 
typically cause the LF to appear to evolve more shal- 
lowly at extreme redshifts. However, the success of our 
selection method in providing a photometric redshift for 
183348 of z — 14.32 that was confirmed by Carniani et 
al. (submitted) provides some evidence that our highest 
redshift candidates could bear out. 

Since the remaining candidates at z > 13.5 have some 
imperfection in their cases, as discussed in 8 3.3, and 
to illustrate the relative impact of the highest-redshift 
objects on our inferences, we consider the impact on the 
LF estimate if we were to ignore the z > 14 objects. 
Removing these objects makes the inferred evolution of 
the LF notably steeper, which we show through the UV 
luminosity density evolution in Figure 17 where the light 
jade region and gray line report the marginalized 16-84% 
credibility interval and median puy, respectively. Since 
the fiducial model assumes an evolution oi (z) that has 
a log-linear dependence on redshift, the puy inferred by 
the model beyond the redshift of our observed sample 
can in principle be artificially inflated by the inferred 
trend at z ~ 12—14. Instead, when removing the z > 14 
objects, we explore a more rapid decline given by 


log oy (z) = log($.,0) x exp [(z — 2o)/h;] - (12) 


This model enables a log-exponential drop in the galaxy 
abundance. Indeed, without the z > 14 objects the in- 
ferred pyy would drop much more rapidly than in the 
fiducial model based on the Main Sample. For reference, 
by z ~ 16 the difference between the two inferences is 
more than an order of magnitude. Of course, given the 
small number statistics, we are also sensitive to the im- 
pact of a single false negative. If any of the remaining 
Auxiliary Sample objects in § 3.3 were to prove out, the 
LF would surely rise. 


5.5. Comparison with Halo Abundance and 
Large-Scale Structure 


The large-scale structure of the Universe is expected to 
present a substantial cosmic variance uncertainty given 
the small size of this field. High-redshift galaxies likely 
live in rare halos of high mass for their epoch, leading 
to a large clustering bias and substantial number den- 
sity fluctuations. To investigate this, we utilize the halo 
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catalog from a cold dark matter simulation performed 
by the Abacus N-body code as part of the AbacusSum- 
mit suite (Garrison et al. 2021; Maksimova et al. 2021). 
This simulation used 6144? particles in a 3005! Mpc 
box, resulting in a particle mass of 1.5 x 10° Mo, and 
a force softening of 21 comoving kpc. Halos were found 
using the CompaSO algorithm (Hadzhiyska et al. 2021). 
While this simulation has high accuracy, we caution that 
the measurement of halo mass always depends on the 
halo-finding algorithm; we focus here on the relative 
trends across redshift and on the clustering. 

In Figure 18, we compare our LF measurements to 
the cumulative halo mass function as a function of red- 
shift. One sees that if the shallow LF is correct, then 
matching the abundance of these galaxies to the abun- 
dance of the most massive halos would require a strongly 
evolving halo mass. On the other hand, if one were to 
discard the objects at z > 14, then the result is more 
similar to the abundance of a constant mass, roughly 
of 101? Ma. Of course, the galaxies may live in less 
massive halos, with a scatter between luminosity and 
mass (e.g., Shen et al. 2023b; Sun et al. 2023); indeed, 
some scatter is inevitable (Pan & Kravtsov 2023). In 
what follows, we therefore consider the properties of ha- 
los with virial masses of 10?" Mọ, about 340 particles, 
which has comparable abundance to our galaxy sample 
at z ~ 12 — 14. 

We then calculate the variation within the simulation 
of regions similar in size to the JOF. We use pencil- 
shaped regions of 6^! Mpc square, roughly 3' at z ~ 12, 
with a depth appropriate to Az — 1. We find that at 
z — 12 (11.5-12.5), there are an average of 8.3 halos 
above our mass threshold in a region, but with a stan- 
dard deviation of 5.6. At z — 13, this abundance drops 
to 2.3 + 2.3; at z = 14, the abundance drops further 
to 0.7167. The distribution of halo number counts be- 
comes noticeably skewed, and by z — 14 we find that 
6% of regions have > 3 halos. Hence, we find that 
unless the host halos are much less massive (and their 
luminosity much more variable), the large-scale struc- 
ture contributes an error at least as large as the Poisson 
error. We caution that this uncertainty could impact 
the observed rate of decline of the UV luminosity den- 
sity, given our area, and motivates further studies over 
larger fields. However, to combat other systematics such 
studies should also leverage the depth and filter cover- 
age comparable to that afforded by the JOF, which is 
challenging given the necessary exposure time. 

Finally, we note that we have neglected the effect of 
magnification by gravitational lensing in our inference of 
Muy. While none of our candidates show obvious lens 
morphology, the high-luminosity tail of the high-redshift 
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Figure 16. UV luminosity at z ~ 12 inferred from the 
JADES Origins Field (JOF). Using the method described in 
85.1, we compute the marginalized constraints on the UV 
luminosity function inferred from galaxies discovered in the 
JOF with photometric redshift distributions that overlap the 
redshift range 11.5 « z « 13.5. We account for photomet- 
ric scatter, the photometric redshift distribution of each ob- 
ject, the selection completeness for each object, and poten- 
tial contamination from proximate redshifts. The 16976-8496 
marginal constraints on the abundance uv as a function 
of absolute UV magnitude Myy are shown as a jade-shaded 
area and the median $vv(Muv) is shown as a white line. 
For comparison, we also compute step-wise luminosity func- 
tion constraints as described in 85.2 at z ~ 12 (solid black 
points) and at z ~ 14 (open black circles). These step-wise 
estimates agree with the inferred ovv, but the continuous 
constraints represent our results for the UV LF. We also 
show a variety of constraints from the literature at compa- 
rable redshifts (colored points), and note that none of these 
data were used to aid our inference of the UV LF. 


luminosity function will likely be enhanced by lensing 
(e.g., Wyithe et al. 2011; Mason et al. 2015; Ferrami & 
Wyithe 2023), which might affect interpretations of the 
luminosity function in the context of theories of galaxy 
formation. 


6. PHYSICAL PROPERTIES OF THE 
HIGH-REDSHIFT POPULATION 


Beyond the abundance and UV luminosity of these 
z = 12 galaxies, the physical properties of the galax- 
ies are of particular interest for understanding the pro- 
cess of galaxy formation at the earliest epochs. With 
the high-quality space-based optical-infrared photome- 
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Figure 17. Evolution of the UV luminosity density ouv(Muv < —17) with redshift derived from the JOF sample. Shown are 
literature values for puv(z) measured (left panel) or extrapolated (right panel) to Muy < —17. In both panels, the shaded 
jade region shows the 16% and 84% marginal constraints on the luminosity density computed from the posterior samples of the 
evolving luminosity function inference, as well as the median luminosity density with redshift (white line). These constraints 
model a linear evolution in log,, x and include a permissive prior on the faint-end slope a. Overall, our constraints agree well 
with prior literature results even as our inference is completely independent. ‘The dark green lines extending to z ~ 8 show 
the low-redshift extrapolation of the inferred puv(z) evolution, while the shaded region indicates the redshift range where our 
detection and selection completeness is non-negligible. We also indicate an approximate cosmic star formation rate density 
(right axis; Mo yr ` Mpc ?) using the conversion syy = 1.15 x 10 ^5 Mo yr ^! erg! s Hz, and show the Madau & Dickinson 
(2014) model (left panel; dotted line). For comparison, inn the left panel, we show the corresponding constraint if the JOF 
high-redshift galaxies and candidates at z > 14 are excluded and log, ©, is fit with an exponential evolution. In this case, we 
would infer the light jade region (1696-8476 marginal constraint) with gray line (median). 


try available in the JOF, physical properties of the high- 
redshift galaxy stellar populations can be inferred. 


6.1. Rest-frame UV Magnitude and Spectral Slope 


Given the dramatic distances to these objects, the 
photometry obtained in the JOF primarily probes only 
their rest-frame UV spectra. Using common-PSF im- 
ages and aperture-corrected Kron photometry as a proxy 
for the total fluxes, we can fit the rest-frame UV pho- 
tometry with a power law f, x A?*? and jointly con- 
strain Myy and 8 given the object redshifts. Figure 19 
shows the posterior distribution of Myy and 8 for the 
candidate galaxies in our Main Sample at z > 11.5. The 
posterior mean and standard deviation for each param- 
eter are reported in Table 8, and for convenience we also 
report Myy in Table 2. The maximum likelihood val- 
ues for the rest-frame spectral slope are —2 = B = —3. 
These values are comparable to the rest-frame spectral 
properties of high-redshift photometric samples (e.g., 
Cullen et al. 2023a; Topping et al. 2023), although not 
quite blue enough to suggest completely dust-free ob- 
jects (e.g., Cullen et al. 2023b). 


6.2. Morphology and Size 


As expected, these galaxies show small angular sizes. 
As described in 2.3.1, we fit single Sérsic profiles to the 
individual exposures in the F200W and F277W filters, 
reporting the half-light radii in Table 2. The posterior 
distributions are often non-Gaussian and asymmetric. 
Unsurprisingly, most of the objects are small, with half- 
light radii below 50 mas, excepting the unusual z = 
14.32 galaxy 183348. 

To characterize the limiting angular resolution of our 
images, we have also fit Sérsic profiles to the exposures 
(separated by epoch of observation) in the same bands 
for known brown dwarfs of similar flux levels in the JOF 
and wider GOODS-S areas (Hainline et al. 2023b). As 
in our past work (Robertson et al. 2023), we find that 
brown dwarfs in the JADES Deep imaging are recovered 
with 95% upper limits on sizes of 20 mas in F200W, so 
we regard objects with 95% lower limits above 20mas 
as inconsistent with a point source. As such, candidates 
16699, 160071, and 55733 are resolved, with half-light 
angular sizes up to 50 mas and half-light physical sizes 
of 132, 118, and 142 pc, respectively. The galaxy 183348 
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Table 8. Sample physical properties, assuming best-fit redshift. 


Name NIRCam ID ` Zphot Muv B logig M. [Mo] SFR [Mo yr] 
JADES+53.09731-27.84714 74977 11.53 —17.66 +0.14 —2.09 + 0.28 (08. O47. 
JADES+53.02618-27.88716 16699 11.56 —17.94+0.15 —2.91 + 0.35 DUNT. none 
JADES+53.04017-27.87603 33309 12.1 —17.73 +0.10 —2.46 + 0.24 621% EHS 
JADES+53.03547-27.90037 160071 12.38 —18.16+0.11 —2.43 + 0.27 Fee 0200 
JADES+53.06475-27.89024 13731 12.93 —18.78 +0.04 —2.73 + 0.13 ri nube dete 
JADES+53.02868-27.89301 11457 13.52 —18.55+0.11 —2.46 + 0.30 EHS E Lie 
JADES+53.07557-27.87268 376946 14.38 —18.30+0.22 —2.42+ 0.56 rA cups 0.96 1-23 
JADES+53.08294-27.85563 183348 14.39 —20.93+ 0.04 —2.40 + 0.12 8.86602 CAE 
JADES+53.10762-27.86013 55733 14.63 —18.54+0.13 —2.52 + 0.36 “60 D BTE 


NoTE—The UV absolute magnitude Myy and rest-frame UV slope f are jointly fit to common-PSF Kron photometry for 
each object. We report here the mean and standard deviation of other posterior distributions for each parameter. The 
star formation rates are averaged over the last 10 Myr of the inferred star formation histories. 
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Figure 18. Comparison of the inferred evolution of the JOF 
galaxy number density n(z) and the abundance of dark mat- 
ter halos in cosmological simulations. Shown are the inferred 
number density constraints (dark jade region 16-84%, white 
line 50%) for model with a linear evolution in log,, d, with 
redshift z. The grid of gray lines show the abundance of dark 
matter halos with masses greater than logio M ~ 9.4 — 11 
computed from the AbacusSummit simulation suite (Mak- 
simova et al. 2021). In the inferred JOF n(z), if simply 
matched by abundance the halo mass of the typical galaxy 
would vary by roughly a factor of — 10. If instead we were 
to discard the z > 14 objects and fit an exponential evo- 
lution to log), x, the typical galaxy would mostly track a 
halo mass = log,, M ~ 10 (light jade region). For reference, 
we indicate the extrapolation of the inferred number density 
constraints to lower redshifts with jade lines. 


spectroscopically-confirmed at z = 14.32 by Carniani et 
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Figure 19. Posterior distributions of rest-frame UV ab- 
solute magnitude Myy and spectral slope 6 for candidate 
galaxies in our Main Sample at z > 11.5. Shown as kernel- 
density-estimated contours are the 68% and 95% credibility 
intervals on the joint posterior distributions for each object. 
The maximum likelihood values for the UV spectral slope 
are -22 B = —3. 

The outlier at Muy & —21 is 183348, spectroscopically con- 
firmed at z — 14.32 (Carniani et al., submitted). 


al. (submitted) shows a size of 76 mas, or about 240 
pc. The remaining sources are consistent with a point 
source, though many have non-negligible probability of 
having larger sizes. We note that objects 13731 and 
316946 are both constrained to be very small. In addi- 
tion to the multiband Forcepho fit reported in Table 2, 


independent single-band Forcepho fits to the 13731 infer 
its size be less than 10 and 16 mas (95th percentile) in 
F200W and F277W respectively. While 376946 appears 
unresolved in F200W and F277W, it appears more ex- 
tended in some medium band filters. Regardless, the 
sizes of these objects are small enough that we expect 
their extents do not impact their detection completeness 
(e.g., Figure 13). 

These results are similar to those found in Robertson 
et al. (2023), where 2 of the 4 z > 10 galaxies were re- 
solved. One consequence of being resolved is that the 
light from these galaxies cannot be purely from an ac- 
creting massive black hole (Tacchella et al. 2023). Other 
spectroscopically-confirmed galaxies at z > 12 have had 
size measurements inferred from scene modeling, and 
show sizes of Rı;2 ~ 100 — 300pc (e.g., Wang et al. 
2023b). Collectively, these results indicate that com- 
pact sizes are a common property of many high-redshift 
galaxies and candidates. 


6.3. Star Formation Rate Histories 


To perform detailed modeling of the SEDs in terms of 
stellar populations, we use the Prospector code (John- 
son et al. 2021), following the methods described in Tac- 
chella et al. (2022, 2023). Briefly, we assume a variable 
star-formation history (SFH) with a bursty continuity 
prior, with 8 time bins spanning 0 — 5 Myr, 5 — 10 Myr 
and 6 logarithmically spaced up to z — 25. We allow the 
redshift to vary within the EAZY posterior. We adopt 
a single metallicity for both stars and gas, assuming 
a truncated log-normal centered on log(Z/Ze) = —1.5 
with width of 0.5, minimum of -2.0, and maximum of 
0.0. We model dust attenuation using a two-component 
model with a flexible attenuation curve. For the stel- 
lar population synthesis, we adopt the MIST isochrones 
(Choi et al. 2016) that include effects of stellar rotation 
but not binaries, and assume a Chabrier (2003) initial 
mass function (IMF) between 0.08 and 120 Mo. No 
Lyo emission line is added to the model to account for 
resonant absorption effects, while the IGM absorption 
model (Inoue et al. 2014; Madau 1995) is taken into ac- 
count (normalization is a free parameter). We do not try 
to constrain independently the effects of possible addi- 
tional Lyman-o damping-wing absorption. For consis- 
tency with Figures 2-10, we use the r — 0.1" aperture 
fluxes, but we note that using r — 0.3" aperture fluxes 
provide quantitatively similar results for these compact 
objects. We put an error floor of 596 on the photometry. 
The rest of the nebular emission (emission lines and con- 
tinuum) is self-consistently modeled (Byler et al. 2017) 
with two parameters, the gas-phase metallicity (tied to 
the stellar metallicity), and the ionization parameter 
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(uniform prior in —4 < log(U) < —1). By combining 
these inferred stellar population properties with the size 
measurements from ForcePho, we can additionally infer 
the stellar mass and star formation rate surface densities 
of the candidate galaxies. 

Figure 20 shows the resulting star formation rate his- 
tories (SFHs) of the eight galaxy candidates in our sam- 
ple. The average SFR over the last 10 Myr is also re- 
ported for each candidate galaxy in Table 8. In each 
case, the continuity prior on the star formation history 
was used to inform the point-to-point star formation rate 
variations in the galaxies. For each object, the photome- 
try listed in Tables 3-5 were used, except for the faintest 
object 74977 (f, ~ 2—3nJy) where the lower SNR Kron 
fluxes were used. We find that the typical star formation 
rate of these objects are SFR z 0.1 — 10 M, yr ^! over 
the last t ~ 10 — 30 Myr. The galaxies formed substan- 
tial fractions of their stars in the recent past, and have 
characteristic ages of just a few tens of millions of years. 
A few of the objects (NIRCam IDs 13731, 33309, 55733, 
74977) show features in their SFHs roughly 10 — 20 Myr 
before their observed epoch, with flat or even falling 
SFR thereafter. We speculate that these features may 
reflect “mini-quenching” events where star formation 
shuts down briefly after exhausting or removing fuel 
(Looser et al. 2023). For the other objects, the SFHs 
appear to increase to the epoch of observation, sug- 
gesting some upswing in the star formation rate and 
luminosities of these objects. In two cases (NIRCam ID 
74977 and 183348) the objects show evidence of com- 
parable or higher star formation rates 100 Myr before 
the observed epoch. For 74977, this early star forma- 
tion would correspond to z ~ 14.2. For 183348, the 
early star formation would potentially start at z ~ 20. 
The uncertainties on the SFH are large, and we cannot 
constrain well the star formation rate before z ~ 15 for 
most objects. Given the physical sizes of the objects of 
Ry /2 ~ 50 — 200pc inferred from the ForcePho analysis, 
the star formation rate surface densities of these objects 
are Xspg ~ 10—100 Mo yr! kpc?. Both the SFR and 
SFR surface densities are comparable to those found by 
Robertson et al. (2023) for spectroscopically-confirmed 
galaxies at z ~ 12 — 13, and consistent with being from 
the same galaxy population. 

The above analysis assumes no luminous contribution 
from an active galactic nucleus. Of course, some of these 
galaxies may possibly host luminous AGN, as have been 
found or suspected in some other high-redshift galaxies 
(e.g., Goulding et al. 2023; Übler et al. 2023; Kokorev 
et al. 2023; Maiolino et al. 2023a,b). AGN emission 
would decrease the inferred stellar emission and require 
a re-assessment of the star formation histories and stellar 


28 


10? 10? 
10! 10! 
LC 10° = 10° 
S >, E 
= 10-1 = = 10-1 
© 102 G 8 10* 
am Cc am 
LL LL. LL 
U) 1073 dp) 0) 10-3 
10*t JADES 53.09731-27.84714 JADES 53.02618-27.88716 10“ E JADES 53.04017-27.87603 
ID 74977; z,=11.5 ID 16699; z, — 11.6 ID 33309; z, — 12.1 
E - -5 
10 109 10! 10? 1 10° 10! 10? i 109 10! 10? 
Time [Myr] Time [Myr] Time [Myr] 
10? 10? 
10! 10! 
LC 109 LC = 10° 
E 5, S 
= 10-1 = = 10-1 
Zei E Š 102 
am cc tc 
LL. LL. LL 
N 10% N) Q) 10-3 
10*t  JADES 53.03547-27.90037 JADES 53.06475-27.89024 10*t JADES 53.02868-27.89301 
ID 160071; z, = 12.4 ID 13731; z, = 12.9 ID 11457; z, — 13.5 
-5 - -5 
i 109 10! 10? 19 109 10! 10? 19 109 10! 10? 
Time [Myr] Time [Myr] Time [Myr] 
10? 10? 
10! 10! 
m 109 = c 109 
E >, E 
= 10-1 = = 1071 
8 107 8 8 102 
am (c tr 
LL. LL. LL 
N 10° n 10-3 
10"t JADES 53.07557-27.87268. JADES 53.10762-27.86013 10“F  JADES 53.08294-27.85563 
ID 376946; z, = 14.4 | ID 55733; z, = 14.6 ID 183348; z, = 14.4 
5 -5 -5 
10 109 10! 10? " 109 10! 10? 19 109 10! 10? 
Time [Myr] Time [Myr] Time [Myr] 


Figure 20. Star formation histories (SFHs) inferred using the Prospector code (Johnson et al. 2021), assuming a continuity 
prior and following the methods described in Tacchella et al. (2023). The galaxy candidates show star formation rates of 
SFR x 0.1—1 Mo yr! over the last ~ 10 Myr, measured backward from the epoch of observation. Roughly half of the objects 
show increasing star formation histories, while the others indicate a peak or burst in their star formation rates roughly 10 Myr 
before the observation epoch. This feature may indicate an episode of “mini-quenching” (Looser et al. 2023) in these objects. 
Only one galaxy indicates a comparable or higher SFR t ~ 100 Myr before the observation epoch, such that no object indicates 
evidence of substantial star formation before z ~ 15. Each galaxy is labeled by both their [RA,Dec] designation, photometric 


redshift, and internal JADES NIRCam ID. 


masses, and possibly the photometric redshifts. We note 
that the fact that some of these galaxies are angularly 
resolved implies that some of the emission is stellar. 


6.4. Stellar Mass Distributions 


Figure 21 presents the marginal stellar mass distri- 
butions inferred from Prospector fits to the observed 
photometry. The posterior samples of the galaxy prop- 
erties were used to produce marginal distributions of 
the stellar mass, following the procedure described in 
Robertson et al. (2023). In agreement with Robertson 
et al. (2023), we find that the stellar masses of these 
z ~ 12 — 15 galaxies are M, ~ 107 — 10? Ma. Given 
the sizes of R,;5 ~ 50 — 200pc we measure from the 
surface brightness profiles, the stellar mass surface den- 
sities of the objects are then X, ~ 10? — 10* Ms pc ?. 
For a self-gravitating system, the dynamical timescale 
is then comparable to the star formation timescale in- 
ferred in $6.3. Overall, in agreement with our previ- 
ous findings in Robertson et al. (2023), these objects 
are consistent with rapidly star-forming, compact galax- 
ies with formation timescales comparable to a few dy- 
namical times. Using the simple abundance matching 
comparison with dark matter halos discussed in 85.5, 
we note that matching to number densities would place 
these objects in Mj, ~ 101? Mo dark matter halos, with 
M, /My ~ 10! — 10 ?, well above the present-day stel- 
lar mass to halo mass relations (e.g., Wechsler & Tinker 
2018). 


7. DISCUSSION 


The luminosity function evolution remains the best 
current indicator of the connection between galaxies, 
dark matter halos, and cosmic reionization at the high- 
est redshifts (for a review, see Robertson 2022). These 
results from the JADES Origins Field provide some new 
insight into the process of high-redshift galaxy forma- 
tion. 

The JOF provides the best currently available data 
for probing faint galaxies at redshifts z > 12, given its 
depth and filter array. Using an area twice the size of 
the Hubble Ultra Deep Field, the JOF area reaches a 
deeper limit (30.2 — 30.5AB) and has fourteen JWST 
filters including the ultradeep JADES Program 1210. 
The inclusion of deep F162M provides an essential check 
on the reality of the highest-redshift candidates. 

Of our Main Sample, none of the galaxies are brighter 
than Myy = —18.6, and many have Myy > —18. The 
depth allows us to constrain the UV luminosity func- 
tion to fainter limits at z ~ 14 than previously possible, 
while retaining tighter control of systematics by hav- 
ing additional medium band filters to probe the Lyman 
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Figure 21. Posterior distribution of stellar mass for candi- 
date z > 11.5 galaxies. Shown are the stellar mass distribu- 
tions constructed from posterior samples of the Prospector 
code (Johnson et al. 2021). The objects have inferred stellar 
masses of M, ~ 10’ — 108 Mo, comparable to that inferred 
for the spectroscopically-confirmed z ~ 12 — 13 analyzed by 
Robertson et al. (2023). Each galaxy candidate is labeled 
by its JADES NIRCam ID and photometric redshift, and 
color-coded the same in Figure 20. 


break with more fidelity. Following the stellar popula- 
tion modeling procedure of Tacchella et al. (2023), we 
find that the star formation rate and stellar mass prop- 
erties are comparable to galaxies spectroscopically con- 
firmed at z ~ 12—13 (Robertson et al. 2023; Curtis-Lake 
et al. 2023; Wang et al. 2023b). Using the ForcePho 
forward model for the surface brightness distribution of 
these galaxies, we find that they have compact sizes of 
F1 j © 50 — 200pc, also in agreement with spectroscop- 
ically confirmed galaxies at these redshifts (Robertson 
et al. 2023; Wang et al. 2023b). 

In agreement with previous determinations of UV lu- 
minosity function in extragalactic JWST fields (McLeod 
et al. 2023; Donnan et al. 2023b; Adams et al. 2023b; 
Harikane et al. 2023b,a; Pérez-González et al. 2023a; 
Willott et al. 2023; Finkelstein et al. 2023a), we find 
that the luminosity function of galaxies has smoothly 
declined from z ~ 8, as first established by HST ob- 
servations (e.g., McLure et al. 2013), to z ~ 12. Our 
results for the abundance of galaxies at z ~ 12 are in 
broad agreement with the literature values, as shown in 
Figures 16 and 17. We do note that our inferred UV 
luminosity density at z ~ 14 is lower than that reported 
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by Finkelstein et al. (2023a), but the uncertainties are 
large. 

However, our selection completeness using the JOF 
observations is sensitive to galaxies out to z ~ 20 when 
the Lyman-a break enters F250M. With a suitable re- 
vision to our selection, we would be sensitive to bright 
galaxies at even greater distances. Our work presents 
a new method for modeling the redshift-dependent UV 
luminosity function incorporating both detections and 
non-detections to constrain its evolution over the red- 
shift range z ~ 11 — 20 where our completeness is high. 
From the lack of galaxy candidates at z > 15, we find 
that the decline to z > 14 continues at dlog ¢,/dz ~ 
—0.2 with our nominal Main Sample presented in Ta- 
bles 2-5. We note that uncertainties owing to cosmic 
variance are clearly non-negligible for the JOF, and a 
larger sample of galaxies at z > 11.5 is needed to con- 
firm this decline. Nonetheless, we now know that the 
Myy ~ —21 object NIRCam ID 183348 selected by our 
JOF Medium band photometry to be at a photometric 
redshift of z & 14.4 has been spectroscopically confirmed 
at z — 14.32 by Carniani et al. (submitted). As Fig- 
ure 18 shows, the evolving luminosity density at z > 14 
we infer from 183348 and our photometric candidates, 
while declining, still requires a constant remapping be- 
tween galaxy and halo abundance, with increasing efh- 
ciency in low-mass halos at higher redshifts. This evo- 
lution is in contrast to the possibility that z > 14 galax- 
ies were not abundant, where a rapid drop in the UV 
luminosity density would track more closely the abun- 
dance of Mi ~ 101? Mo halos and the galaxy efficiency 
could stabilize at early times. Given the confirmation of 
183348, we see no evidence for such a stabilization in the 
efficiency of galaxy formation out to z ~ 14 or beyond. 

Lastly, since our results are consistent with prior lit- 
erature results at z ~ 12, theoretical models that match 
those observations also match ours. For instance, the 
feedback-free models of Dekel et al. (2023) and Li et al. 
(2023) agree with our z ~ 12 observations for an effi- 
ciency of Emax © 0.2. Models for the evolving number 
counts of high-redshift galaxies based on dust-free pop- 
ulations (e.g., Ferrara et al. 2023) also predict a star 
formation rate density evolution to z ~ 15 in agreement 
with our inferences, assuming all our candidates are re- 
ally high-redshift sources (Ferrara 2023). 


8. SUMMARY AND CONCLUSIONS 


Using ultra-deep JWST observations of the JADES 
Origins Field (JOF), we search for the most distant 
galaxies in the universe. With fourteen JWST and up 
to nine Hubble Space Telescope filters covering the JOF, 
we can carefully select galaxies at z > 12 by identifying 


dropouts in NIRCam F162M and bluer filters using SED 
template-based photometric redshift fitting. Our find- 
ings include: 


e We select nine galaxy candidates at z ~ 12 — 15 
and no galaxy candidates at z = 15. These objects 
include the most distant candidates detected in 
more than five filters and displaying a dropout in 
more than 10 filters. Our sample selection includes 
a galaxy at z — 14.32 since spectroscopically con- 
firmed. Simulations of our detection and photom- 
etry methods and our prior spectroscopic confir- 
mations of high-redshift JADES sources suggest 
that the other candidates without spectroscopic 
confirmation are robust. Several of our candidates 
have been identified in previous analyses, includ- 
ing Hainline et al. (2023a) and Williams et al. 
(2023). 


These objects show apparent total magnitudes of 
map ~ 29.5 — 30.5 in the rest-frame UV and blue 
rest-UV spectral slopes —2 = 8 = —3. 


Performing detailed structural modeling with For- 
cePho and stellar population inference using 
Prospector, we find that the galaxies have star- 
formation rates of SFR & 0.1 — 10 Mo yr 1, stel- 
lar masses of M, ~ 107 — 10? Ms, sizes of R ~ 
50 — 200 pc, and stellar ages of t, ~ 30 — 50 Myr. 
The properties of our low-mass candidates are 
comparable to the properties of z ~ 12 — 13 galax- 
ies with confirmed redshifts, as first identified by 
the JADES collaboration. 


e We develop à new forward modeling method to 
infer constraints on the evolving UV luminosity 
function without binning in redshift or luminosity 
while marginalizing over the photometric redshift 
posterior distribution of candidates in our sample. 
This method allows for an accounting of poten- 
tial contamination by adjacent redshifts and in- 
cludes the impact of non-detections on the inferred 
galaxy luminosity function evolution. 


With the population of z > 12 galaxy candidates 
newly discovered in JOF, we provide an inference 
on the z ~ 15 luminosity function and a refined 
measure of the luminosity function at z ~ 12 in 
agreement with literature values. At z ~ 15, we 
infer a continued decline from z ~ 12. Over the 
redshift range z ~ 12 — 14, where we have de- 
tected galaxies, we infer a factor of 2.5 decline 
in the luminosity function normalization $, and 
a corresponding decline in the luminosity density 


puv. We note that cosmic variance uncertainties 
for the high-redshift JOF sample are not negligi- 
ble, and this decline should be confirmed with a 
larger sample over a wider area. 


This demonstrates the immediate impact new JWST 
observations can have on our knowledge of the distant 
universe. With high-redshift galaxy populations now 
established fewer than 300 million years after the Big 
Bang, we have extended our reach into the cosmic past 
by 40% during the first eighteen months of JWST oper- 
ations. 
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